Multiple physical elements to determine the gravitational-wave 
signatures of core-collapse supernovae 



D4 



Kei Kotake''''^ 

"Division of Theoretical Astronomy, National Astronomical Obserx'atory of Japan, 2-21-1, Osawa, Mitaka, Tokyo, 

181-8588, Japan 

Center for Computational Astrophysics, National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 



Abstract 

We review recent progress in the theoretical predictions of gravitational waves (GWs) of core- 
collapse supernovae. Following a brief summary of the methods in the numerical modeling, 
we summarize multiple physical elements that determine the GW signatures which have been 
considered to be important in extracting the information of the long-veiled explosion mechanism 
from the observation of the GWs. We conclude with a summary of the most urgent tasks to make 
the dream come true. 
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1. Introduction 

Massive stars in the approximate mass range of 8 to 100 solar masses (Mq) end their lives 
as core-collapse supernovae. They have long attracted the attention of astrophysicists because 
they have many aspects playing important roles in astrophysics. They are the mother of neutron 
stars and black holes; they play an important role for acceleration of cosmic rays; they influence 
galactic dynamics triggering further star formation; they are gigantic emitters of neutrinos and 
gravitational waves. They are also a major site for nucleosynthesis, so, naturally, any attempt 
to address human origins may need to start with an understanding of core-collapse supernovae 
(SNe). 

Current estimates of core-collapse SN rates in our Galaxy predict one event every ~ 40 ± 10 
year f!]]. When a massive star undergoes a core-collapse SN in our Galactic center, copious num- 
bers of neutrinos are produced, some of which may be detected on the earth. Such "supernova 
neutrinos" will carry valuable information from deep inside the core. In fact, the detection of 
neutrinos from SN1987A (albeit in the Large Magellanic Cloud) opened up the neutrino astron- 
omy, which is an alternative to conventional astronomy by electromagnetic waves Even 
though there were just two dozen neutrino events from SN1987A, these events have been stud- 
ied extensively (yielding ~ 500 papers) and have allowed us to have a confidence that our basic 
picture of the supernova physics is correct (e.g., see jst] for a recent review). If a supernova 
occurs in our Galactic center (~ 10 kpc), about 10,000 Ve events are estimated to be detected 
by Super- Kamiokande (SK) |@] (e.g., jTiHtl). Those successful neutrino detections are important 
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not only to study the supernova physics but also to unveil the nature of neutrinos itself such as 
the neutrino oscillation parameters and the mass hierarchy (e.g., |§| for a recent review). 

Core-collapse SNe are now about to start even another astronomy, Gravitational-Wave As- 
tronomy. Currently long-baseline laser interferometers such as LIGO (USA)IToll, VIRGO (Italyfl 
GEO600 (Germany jl, and TAMA300 (Japan) [[Till are currently operational and preparing for the 
first observation (see, e.g., lll2ll for a recent review), by which the prediction by Einstein's the- 
ory of General Relativity (GR) can be confirmed. These instruments are being updated to their 
Advanced status, and may start taking data, possibly detecting GWs for the first time, as soon as 
2015 (see lll3ll for a recent review). In fact, AdvancedLlGOfVlRGO, which is an upgrade of the 
initial LIGO and VIGRO, are expected to be completed by 2015 and will increase the observ- 
able detection volume by a factor of ~ 1000 il4ll . The Large-scale Cryogenic Gravitational-wave 
Telescope (LCGT ifTsIl ) in Japan was funded in late 2010, which is being built under the Kamioka 
mine and is expected to take its first data in 2016. At such a high level of precision, those GW 
detectors are sensitive to many different sources, including chirp, ring-down, and merger phases 
of black-hole and neutron star binaries (e.g., l3i 18]), neutron star normal mode oscillations 
(e.g., 19 ]), rotating neutron star mountains (e.g., j20tl), and core-collapse supernova explosions 



(e.g., ll2lLl22il23ll for recent reviews), on the final of which we focus in this article. 



According to the Einstein's theory of general relativity (e.g., 112411 '). no GWs can be emit- 
ted if gravitational collapse of the supernova core proceeds perfectly spherically symmetric. To 
produce GWs, the gravitational collapse should proceed aspherically and dynamically. Observa- 
tional evidence gathered over the last few decades has pointed towards core-collapse SNe indeed 
being generally aspherical (e.g., ||25, 2^ 27, 28] and references therein). The most unequivocal 
example is SN1987A. The HST images of SN1987A are directly showing that the expanding 
envelope is elUptical with the long axis aligned with the rotation axis inferred from the ring 
( ll26ll . see however ll29ll for a recent counter argument). The aspect ratio and position angle of the 
symmetry axis are consistent with those predicted earlier from the observations of speckle and 
linear polarization. What is more, the linear polarization became greater as time passed (e.g., 
lEs . 30 . 31]), a fact which has been used to argue that the central engine of the explosion is re- 
sponsible for the non-sphericity (e.g., ]l32, 33]). From a theoretical point of view, clarifying what 



makes the dynamics of the core deviate from spherical symmetry is essential in understanding the 
GW emission mechanism. Here it is worth mentioning that GWs are primary observables, which 
imprint a live information of the central engine, because they carry the information directly to 
us without being affected in propagating from the stellar center to the earth. On the other hand, 
SN neutrinos, which are affected by the well-known Mikheyev-Smirnov-Wolfenstein (MSW) 
effect (e.g., 13411 ) in propagating the stellar envelope, may be rather indirect in extracting the 
information of the central engine compared to GWs|j 

The breaking of the sphericity in the supernova engine has been also considered as a most 
important ingredient to understand the explosion mechanism, for which supernova theorists have 
been continuously keeping their efforts for the past ~ 40 years. Albeit being a topic of hot debate 
yet, current multi-dimensional (multi-D) simulations based on refined numerical models show 
several promising scenarios. Among the candidates are the neutrino heating mechanism aided 
by convection and hydrodynamic instabilities of the supernova shock (e.g., ll36ll for a review). 



' http://www.ego-gw.it/ 
'http://geo600.aei.mpg.de/ 

^However they can be a useful tool like a tomography to monitor the density profile in the stellar mantle, which is 
expected to tell us the evolution of the supernova shock (e.g., for a recent review). 
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the acoustic mechanism OTIl . or the magnetohydrodynamic (MHD) mechanism (e.g., 112 ih see 
references therein). 

Putting things together, the multi-dimensionahty determines the explosion mechanism, in 
turn we may extract the information that traces the muhi-dimensionaHty by the SN messengers, 
which would be only possible by a careful analysis primary on GWs as well as neutrinos, and 
photons. In this review article, we hope to bring together findings from recent multi-D supernova 
simulations and the obtained predictions of the GWs so far (for GWs from other high-energy 
astrophysical sources such as magnetars, gamma-ray bursts, and coalescing binaries, see 138 , 
39, ^ 41] for recent reviews). Before moving on to the next sections, we first have to draw 



a caution that the current generation of numerical results that we report in this article and the 
associated predictions of the GW emission should depend on the next-generation calculations 
by which more sophistication can be made not only in determining the efficiency of neutrino- 
matter coupling (the so-called neutrino transport problem), but also in the treatment of general 
relativity. Therefore we provide here only a snapshot of the moving (long-run) documentary 
film that records our endeavours for making our dream of "GW astronomy of core-collapse 
supernovae" come true. 



2. On core-collapse supernova simulations 

In order to extract a true astrophysical GW signal from the confusing detector noises, we 
need not only require sensitive detectors, but also extensive data analysis of the detector output 
based on reliable theoretical estimates for the GW signals. In most cases, a detailed numerical 
modeling is required to make the precise prediction of the dynamics for the GW sources. As will 
be summarized briefly in this section, SN modelers have been keeping their efforts to this end for 
these -1-45 years. 

The abstract in the paper by Colgate and White in 1966 I42il who reported the first SN simu- 
lation finished with the sentence, "77ie energy release (at the moment of explosion) corresponds 
to the change in gravitational potential of the unstable imploding core; the transfer of energy 
takes place by the emission and deposition of neutrinos'^ As well-known, this is the essence of 
the so-called neutrino-heating mechanism of core-collapse supernovae. The mechanism was first 
proposed to occur in a prompt manner, but was later reinforced by Bethe and Wilson ll43il44ll45ll 
to take a currently prevailing delayed form, in which a stalled bounce shock is revived by neu- 
trino energy deposition to trigger explosions in several hundred milliseconds after bounce. The 
mechanism has been the working hypothesis of supernova theorists for these ~45 years (see col- 
lective references in Janka et al. (2007) Isill that was presented in the Bethe Centennial Volume 
in Physics Report, and also a summary given in section 2 of Nordhaus et al. (2010) ll46ll ). 

However, one important lesson we have learned from a series of the most up-to-date sim- 



ulations 143 |48l |49l |50tl 

which implemented the best input physics and numerics to date, is 



that the mechanism fails to blow up canonical massive star^in spherical symmetric (ID) sim 
ulations. Pushed by supernova observations of the blast morphology [e.g., |25l 12811 mentioned 
above, it is now almost certain that the breaking of the spherical symmetry is the key to solve 
the supernova problem. So far a number of multi-D hydrodynamic simulations have been re- 
ported, which demonstrated that hydrodynamic motions associated with convective overturn 



''l boldly added several words in (...) for making the meaning of the sentence complete, 
'except for 8 - 9Mq stars (e.g., (sTj]). 
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5a\) as well as the Standing-Accretion-Shock-Instability (SASI, e.g., 
63l 164 l65l l66l l67l 16811 see references therein) can help the onset of 



the neutrino-driven explosion. 

To test the neutrino heating mechanism in the multi-D context, it is of crucial importance 
to solve accurately the neutrino-matter coupling in spatially non-uniform hydrodynamic envi- 
ronments. For the purpose, one ultimately needs to solve the six-dimensional (6D) neutrino 
radiation transport problem (three in space, three in the momentum space of neutrinos), thus 
the supernova simulation has been counted as one of the most challenging tasks in numerical 
astrophysics. In the final sentence of the last paragraph, we wrote "demonstrated" because the 
neutrino heating was given by hand as an input parameter in most of the simulations cited above 
(see, however 1551 15611 ). The neutrino heating proceeds dominantly via the charged current in- 
teractions (ve -H n e -H p, Ve -H p ^ e^ H- n) in the gain region. The neutrino heating rate 
in the gain region can be roughly expressed as 16911 oc Ly{ijyy^ /r^ where Ly is the neutrino 
luminosity emitted from the surface of neutrino sphere and it determines the amplitude of the 
neutrino heating as well as coohn g, a nd r and {fiy} is the distance from the stellar center and 
the flux factoid, respectively (e.g., 1691). For example, L,, is treated as an input parameter in the 
so-called 'light-bulb" approach (e. g., 5411 ). This is one of the most prevailing approximations in 
recent 3D simulations 1661 67 . 70l 46 1 because it is handy to study multi-D effects on the neu- 
trino heating mechanism (albeit on the qualitative grounds). To go beyond the light-bulb scheme, 
Ly should be determined in a self-consistent manner For the purpose, one needs to tackle with 
neutrino transport problem, only by which energy as well as angle dependence of the neutrino 
distribution function can be determined without any assumptions. Since the focus of this review 
is on GWs, a detailed discussion of various approximations and numerical techniques taken in 
the recent radiation-hydrodynamic SN simulations cannot be provided. Table 1 is not intended 
as a comprehensive compilation, but we just want to summarize milestones that have recently 
reported exploding models so far The table will be useful in later sections when we discuss the 
GW signatures obtained in some of these simulations. 

In Table l,the first column ("Progenitor") shows the progenitor model employed in each sim- 
ulation. The abbreviation of "NH", "WHW", and "WW" means Nomoto & Hashimoto (1988) 
III, Woosley, Heger, and Weaver (2002) I?!, and Woosley and Weaver (1995) H- The sec- 
ond column shows SN groups with the published or submitted year of the corresponding work. 
"MPA" stands for the core-collapse SN group in Max Planck Institute for Astrophysics led by 
H.T. Janka and E. Miiller. "Princeton-H" stands for the group chiefly consisting of the staffs in 
the Princeton University (A. Burrows), Caltech (C.D.Ott), Hebrew University (E. Livne) and 
their collaborators. The SN group in the Basel university is led by M. Liebendorfer and FK. 
Thielemann. "OakRidge+" stands for the SN group mainly consisting of the Florida Atlantic 
University (S. Bruenn) and the Oak Ridge National Laboratory (A. Mezzacappa, O.E.B. Messer) 
and their collaborators. NAOJ+ is the SN group chiefly consisting of the staffs in the National 
Astronomical Observatory of Japan (myself, T. Takiwaki), Kyoto university (Y. Suwa), Waseda 
University (S. Yamada), and their collaborators. The third column represents the mechanism of 
explosions which are basically categorized into two (to date), namely by the neutrino-heating 
mechanism (indicated by "y-driven") or by the acoustic mechanism ("Acoustic"). "Dim." in the 
fourth column is the fluid space dimensions which is one-, two-, or three-dimension (1,2, 3D). 
The abbreviation "N" stands for 'Newtonian,' while "PN" — for 'Post-Newtonian' — stands for 



*This quantity represents the degree of anisotropy in neutrino emission; (/ly) 
(/iy) = 1 in the free-streaming limit (r — * oo) (e.g., Janlca (2001) |6^ ). 
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~ 0.25 near at the neutrino sphere. 



some attempt at inclusion of general relativistic effects, and "GR" denotes full relativity, fexp in 
the fifth column indicates an approximate typical timescale when the explosion initiates and ^exp 
represents the explosion energy normalized by Bethe(=10^' erg) given at the postbounce time of 
fpb, both of which are attempted to be sought in literature^. In the final column of "v transport", 
"Dim" represents v momentum dimensions and the treatment of the velocity dependent term in 
the transport equations is symbolized by 0(v/c). The definition of the "RBR", "IDSA", and 
"MGFLD" will be given soon in the following. 

Due to the page limit of this article, we have to start the story only after 2006 (see, e.g., 
Janka et al. (2007) ll36ll for a grand review, and also Cardall (2005) 1 83 1 for a similar table before 
2006). The first news of the exploding model was reported by the MPA group. By using their 
MuDBaTH code which includes one of the best available neutrino transfer approximations, they 
reported ID and 2D explosions for the 8.8 Mq star by NH88 whose progenitor has a very tenu- 
ous outer envelope with steep density gradient, which is a characteristic property of AGB stars. 
Also in ID, the Basel+ group reported explosions for 10 and 15 Mq progenitors of WHW02 
triggered by the hypothesised first-order QCD phase transition in the protoneutron star (PNS). 
To date, these two are the only modern numerical results where the neutrino-driven mechanism 
succeeded in ID. In the 2D MPA simulations, they obtained explosions for a non-rotating 11.2 
Mq progenitor of WHW02 01 , and then for a 15Mo progenitor iH of WW95 with a relatively 
rapid rotation imposecfl They newly brought in the so-called "ray-by-ray" approach (indicated 
by "RBR" in the table), in which the neutrino transport is solved along a given radial direction as- 
suming that the hydrodynamic medium for the direction is spherically symmetric. This method, 
which reduces the 2D problem partly to ID, fits well with their original ID Boltzmann solver 
1 47 E For 2D hydrodynamic simulations with the ray -by-ray transport, one needs to solve the 4D 
radiation transport problem (two in space and two in the neutrino momentum space). Regarding 
the explosion energies obtained in the MPA simulations, their values at their final simulation 
time are typically underpowered by one or two orders of magnitudes to explain the canonical 
supernova kinetic energy (~ 10^' erg). But the explosion energies presented in their figures are 
still growing with time, and they could be as high as 1 B if they were able to follow a much 
longer evolution as discussed in 17711 . 

More recently, fully 2D multi-angle Boltzmann transport simulations become practicable 
by the Princeton+ group |8^ 87 1. In this case, one needs to handle the 5D problem for 2D 
simulations (two in space, and three in the neutrino momentum space). However this scheme 
is very computationally expensive currently to perform long-term supernova simulations. In 
fact, the most recent 2D work by IItIi succeeded in following the dynamics until ~ 400 ms 
after bounce for a non-rotating and a rapidly rotating 20 Mq model of WHW02, but explosions 
seemingly have not been obtained in such an earUer phase either by the neutrino-heating or the 
acoustic mechanism. 

In the table, "MGFLD" stands for the Multi-Group Flux-Limited Diffiision scheme which 
eliminates the angular dependence of the neutrino distribution function (see, e.g., Bruenn (1985) 
11881] for more details). For 2D simulations, one needs to solve the 3D problem, namely two in 
space, and one in the neutrino momentum space. By implementing the MGFLD algorithm to the 



^but if we cannot find them, we remain tliem as blank 

^by comparing the precollapse angular velocity to the one predicted in a recent stellar evolution calculation f84ll . 

'Note that the ray-by-ray approach has an advantage compared to other approximation schemes, such that it can fully 
take into account the available neutrino reactions (e.g., fsSl for references therein) and also give us the most accurate 
solution for a given angular direction. 
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Progenitor 


Group 


Mechanism 


Dim. 


^exp 


£exp(B) 


V transport 




(Year) 




(Hydro) 


(ms) 


@fpb (ms) 


(Dim, Oivic)) 




MPA|;51,74j 


v-driven 


1D/2D 


-200 


0.1 


Boltzmann 


8.8 Mo 


(2006) 




(PN) 




(-800) 


2, OivIc) 


(NH88[7]J) 


Princeton+ 


v-driven 


2D 


^125 


0.1 


MGFLD 




|75J(2006) 




(N) 




- 


1,(N) 


10 Mo 


Basel[76] 


v+(QCD 


ID 


255 


0.44 


Boltzmann 


(WHW02[72]) 


(2009) 


transition) 


(GR) 




(350) 


2, (GR) 


11 Mo 


Princeton+ 


Acoustic 


2D 


^550 


-0.1* 


MGFLD 


(WW95[73J) 


[75J(2006) 




(N) 




(1000) 


1,(N) 




MPAI77J 


v-driven 


2D 


-100 


- 0.005 


"RBR" Boltz- 


11.2 Mo 


(2006) 




(PN) 




(-220) 


mann, 2, 0{v/c) 


(WHW02I72J) 


Princeton+ 


Acoustic 


2D 


>1100 


-0.1* 


MGFLD 




[78J (2007) 




(N) 




(1000) 


1,(N) 




NAOJ+ 


v-driven 


3D 


-100 


0.01 


IDSA 




|79j(2011) 




(N) 




(300) 


1,(N) 


12 Mo 


Oak Ridge + 


v-driven 


2D 


-300 


0.3 


"RBR" MGFLD 


(WHW02[72j) 


[80] (2009) 




(PN) 




(1000) 


hO{v/c) 


13 Mo 


Princeton+ 


Acoustic 


2D 


>1100 


-0.3* 


MGFLD 


(WHW02|72j) 


[78j(2007) 




(N) 




(1400) 


1,(N) 


(NH88I71J) 


NAOJ+ 


v-driven 


2D 


-200 


0.1 


IDSA 




[81J(2010) 




(N) 




(500) 


1,(N) 


15 Mo 


MPA[82J 


v-driven 


2D 


-600 


0.025 


Boltzmann 


(WW95[73J) 


(2009) 




(PN) 




(-700) 


2,0(v/c) 


(WHW02I72J) 


Princeton+ 


Acoustic 


2D 






MGFLD 




[78J 




(N) 




(-) 


1,(N) 




OakRidge+ 


v-driven 


2D 


-300 


-0.3 


"RBR" MGFLD 




[801(2009) 




(PN) 




(600) 


l,0{v/c) 


20 Mo 


Princeton+ 


Acoustic 


2D 


>1200 


-0.7* 


MGFLD 


(WHW02|72j) 


[78j(2007) 




(N) 




(1400) 


1,(N) 


25 Mo 


Princeton+ 


Acoustic 


2D 


>1200 




MGFLD 


(WHW02[72j) 


[78J(2007) 




(N) 




(-) 


1,(N) 




Oak Ridge + 


v-driven 


2D 


-300 


-0.7 


"RBR" MGFLD 




[80] (2009) 




(PN) 




(1200) 


l,0(v/c) 



Table 1 : Selected lists of recent neutrino-radiation hydrodynamic milestones reported by many SN groups around the 
world ("Group"), which obtained explosions by the neutrino-heating mechanism (indicated by "y-driven") or the acoustic 
mechanism ("Acoustic") (See text for more details). 
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CHIMERA code in a ray-by-ray fashion (e.g., MW), Bmenn et al. (2009) 18911 obtained neutrino- 
driven explosions for non-rotating progenitors in a relatively wide range in 12, 15, 20, 25 Mq of 
WHW02 (see table). These models tend to start exploding at around 300 m after bounce, and the 
explosion ener gy f or the longest running model of the 25 progenitor is reaching to 1 B at 1 .2 
s after bounce ISOll . 

On the other hand, the 2D MGFLD simulations implemented in the VULCAN code JtsIi 
obtained explosions for a variety of progenitors of 11, 11.2, 13, 15, 20, and 25 Mq not by the 
neutrino-heating mechanism but by the acoustic mechanisrr0- The acoustic mechanism relies 
on the revival of the stalled bounce shock by the energy deposition via the acoustic waves that 
the oscillating protoneutron stars (PNSs) would emit in a much delayed phase (~ 1 second) 
compared to the conventional neutrino-heating mechanism (~ 300 - 600 milliseconds). If the 
core pulsation energy given in Burrows et al. (2007) 17811 could be used to measure the explosion 
energy in the acoustic mechanism, they reach to 1 B after 1000 ms after bounc4l3- The additional 
energy input from acoustic waves is very appealing, but it may remain a matter of vivid debate 
and has yet to be confirmed by other groups. 

By performing 2D simulations in which the spectral neutrino transport was solved by the 
isotropic diffusion source approximation (IDSA) scheme ll90ll . the NAOJ-H group reported ex- 
plosions for a non-rotating and rapidly rotating 13 Mq progenitor of NH88. They pointed out 
that a stronger explosion is obtained for the rotating model comparing to the corresponding non- 
rotating model. The IDSA scheme splits the neutrino distribution into two components (namely 
the streaming and trapped neutrinos), both of which are solved using separate numerical tech- 
niques (see, Liebendorfer et al. (2009) |93l for more details). The approximation level of the 
IDSA scheme is basically the same as the one of the MGFLD. The main advantage of the IDSA 
scheme is that the fluxes in the transparent region can be determined by the non-local distribu- 
tion of sources rather than the gradient of the local intensity like in MGFLD. A drawback in 
the current version of the IDSA scheme is that heavy lepton neutrinos (Vj, i.e., v^, and their 
anti-particles) as well as the energy-coupling weak interactions have yet to be implemented. Ex- 
tending the 2D modules in Suwa et al. (2009) 18111 to 3D, they recently reported explosions in 
the 3D models for an 1 1 .2 Mq progenitor of WHW02 (Takiwaki et al.(201 1) H). By compar- 
ing the convective motions as well as neutrino luminosities and energies between their 2D and 
3D models, they pointed out whether 3D effects would help explosions or not is sensitive to the 
employed numerical resolutions. They argued that next-generation supercomputers are at least 
needed to draw a robust conclusion of the 3D eff'ects. 

Having summarized a status of the current supernova simulations, one might easily see a 
number of issues that remain to be clarified. First of all, the employed progenitors usually rather 
scatter (e.g.. Table 1). Different SN groups seem to have a tendency to employ different pro- 
genitors, providing different results. By climbing over a wall which may have rather separated 
exchanges among the groups, a detailed comparison for a given progenitor needs to be done se- 
riously in the multi-D results (as have been conducted in the Boltzmann ID simulations between 
the MPA, Basel-H, and Oak Ridge-H groups ll9lll ). 

In addition to the importance of 3D modeling as mentioned above, a more complete "realis- 
tic" supernova model should naturally include general relativity (GR) with magnetohydrodynam- 
ics (MHD) with multi-D GR Boltzmann neutrino transport, in which a microphysical treatment 



'"For the 8.8 Mq progenitor, they obtained neutrino-driven explosions (see table 1). 

' ' Due to the ambiguity of the relation between the pulsation energy and the explosion energy, an asterisk is added in 
the table. 



of equation of state (EOS) and nuclear-neutrino interactions are appropriately implemented. Un- 
fortunately none of the currently published SN simulations satisfy the "ultimate" requirement. 
In this sense, all the mentioned studies employ some approximations (albeit with different levels 
of sophistication) for multiple physical ingredients (for example, as listed in Table iF^. 

In the same way, theoretical predictions of the GWs that one can obtain by analyzing the 
currently available numerical results, cannot unambiguously give us the final answer yet. Again 
we hereby note the feature of this article which shows only a snapshot of the moving theoretical 
terrain. Keeping this caveat in mind, it is also true that a number of surprising GW features 
of core-collapse SNe have been reported recently both by the first-principle simulations (e.g., 
in Table 1), and also by idealized simulations in which explosions are parametrically initiated 
mostly by the light-bulb scheme. As will be mentioned in the next section, the latter approach 
is also useful to get a better physical understanding of the GW signatures obtained in the first- 
principle simulation^ We are now ready to move on to focus on the GW signatures from the 
next section. 



3. GW signatures 



We first start to present a short overview concerning possible emission cites of GWs in core- 
collapse SNe so far proposed by a number of extensive studies (section [3.0.1l see also 112 ll |22 . 
2311 for recent reviews). The most up-to-date results regarding the GW signatures emitted near 
bounce in the case of rapidly rotating core-collapse, have been already given in the recent review 
by Ott (2009) (e.g., section 4 in 12211 ). So in this article, we mainly focus on the GWs in the 
postbounce phase in the context of neutrino-driven explosion models and also on the postbounce 
GWs expected in MHD explosions, which will be described separately in sections ITTI and [T2] 



3.0.1. A short overview 

The paper in 1982 by E. Miiller of the MPA entitled as ''Gravitational Radiation from Col- 
lapsing Rotating Stellar Cores" unquestionably opened our eyes to the importance of making 
the GW prediction based on realistic SN numerical modeling (see e.g., section 2 in ll22ll for a 
summary of more earlier work which had mainly focused on the GW emission in very ideal- 
ized systems such as in homogeneous spheroids and ellipsoids0 As one may expect from the 
title of his paper, rapid rotation, if it would exist in the precollapse iron core, leads to signifi- 
cant rotational flattening of the collapsing and bouncing core, which produces a time-dependent 
quadrupole (or higher) GW emission. Following the first study by M iiller , most studies of the 
past thirty have focused on th e so-called bounce signals (e.g., 11981199 . fiMfiMfiolfiolfiol 
lol[io3,[io3,[ioi[Ioi[ii3IIll. 1121) and references therein). 

As summarized by Ott (2009) l22ll . a number of important progresses have been recently 
made to understand features of the bo unce signals by extensive 2D GR studies using the conf orma l 
flat-condition (CFC) approximation iflOSi [l09[ [Tioll and also by fufly GR 3D simulations lfl05 1 
both including realistic EOSs and a deleptonization effect based on ID-Boltzmann simulations 
1 1 1 311 . Due to the page limit of this article, we are only able to touch on them in section l3T2l 



However it is also wortli mentioning that a rash of nice work to this end has been reported recently including new 
schemes towards the GR radiation-hydrodynamic simulations l92Ll93ll94il and new sets of supernova EOSs I95ll96ll97ll . 
In this sense, these two approaches are complimentary in understanding the GW signatures. 
'''Needless to say, this kind of approach is still very important to extract the physics of the GW emission mechanism. 
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For the bounce signals having a strong and characteristic signature, the iron core must ro- 
tate enough rapidly. Although the role of rapid rotation in combination with magnetic fields is 
attracting great attention as an important key in understanding the dynamics of coUapsars and 
magnetars (e.g., section [J!2b . recent stellar ev olutio n calculations predict that such a extreme 
condition can be realized only in a special case 111 14ll^% ~ 1% of massive star population). In ad- 



dition, the precollapse rotation periods are estimated to be larger than ~ 100 sec ll 16ll to explain 
the observed rotation periods of radio pulsars. In such a slowly rotating case, the detection of 
the bounce signals beco mes very hard even by next-generation detectors for a supernova in our 



Galaxy (~ 10 kpc, e.g., 11031110511 ') 



Besides the rapid rotation, anisotropic matter motions associated with convection and anisotropic 
neutrino emission in the postbounce phase are expected to be primary GW sources with compa- 
rable amplitudes to the bounce signals. Thus far, various physical ingredients have been studied 
for producing the post-bo unce asphericit ies, including postshock convection ll 1711, pre -collapse 
density inhomogeneities 1 118 , 1 17 , 1 19l 56ll . moderate rotation of the iron core 1 120ll , non ax- 



isymmetric rotational instabilities 1|12IL 112211 . g-modes 1)12311 and r-modes puls ations III 2411 of 



PNS s, an d more recendy by the Standing-Accretion-Shock InstabiUty (SASI, 01251 Il26lll27 . 

18, I29II V In the mult i-D modeling of stellar evolution, significant progresses have been re- 
cently made (e.g., il30ll and references therein), however the degree of initial inhomogeneities 
in the iron core seems still uncertain. Studies of non-radial instabilities in PNSs and cold NSs 
with their emission processes of GWs have a long history (e.g., section 7 in 12211 for references 
therein). The numerical studies to this end generally treat the microphysics in a very phenomeno- 
logical manner (like by a polytrope EOS), so there may remain a further room for sophistication. 
Among the candidate GW emission mechanisms, we therefore choose to focus on relatively well 
understood parts at first in the next section (section 3.1), which is the GW produced by convection 
and anisotropic neutrino emission. 



3.1. Gravitational waveforms in neutrino-driven explosions 

As mentioned, the first-principle and experimental simulations have a complimentary role to 
understand the GWs in the context of neutrino-driven SN explosions. Since the latter approach 
is equivalently useful, we firstly summarize recent findings obtained by the experimental simu- 
lations in the next sections [3.1.11 and 13 . 1 .21 then proceed to the first-princi ple r esults in section 
13.1.31 For the matter signals, we review the results by Murphy et al. 012911 who conducted 

a systematic study by changi ng ne utrino luminosity and progenitor models. For the neutrino GW 
signals, a piece of our work 112611 is summarized in which a ray-tracing analysis was performed 
to accurately estimate the neutrino GWs. 

3.1.1. Features of Matter GWs in experimental 2D simulations 

Figur e [1] shows samples of gravitational waveforms obtained in 2D models by Murphy et 
al. 1I129II . In their 2D models, the dynamics of a suite of progenitor models (12,15,20, and 40 
Mq) is followed starting from gravitational collapse, through bounce, up to parametric explo- 
sions via the light-bulb scheme. In each panel of Figure [1] the taken progenitor and the input 
(electron-)neutrino luminosity in unit of lO^^erg/s are indicated. For example, top left panel 



'^which experiences the so-called chemically homogeneous evolution IllSll . 

'*Note that "Matter GW" means the GW produced by quadrupole matter deformation and that "neutrino GW" by 
anisotropic neutrino emission. 
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Figure 1 : Sam ples of gravitational waveforms in 2D parametric explosion models (by courtesy of Muiphy and the 
coauthors Il29ll ). Top left panel shows matter GWs obtained for a 15 Mq progenitor model fv^l with an input neutrino 
luminosity of L,,_, = 3.7 X lO'^erg/s. Three typical GW emission cites are clearly indicated, which are "Prompt convection 
(after bounce up to ~ 50 ms after bounce )", "Postshock convection" (then till 550 ms after bounce, unless the word 
is shown in the plot), "Non-linear SASI with SASI plumes"(then till ~ 800 ms after bounce), and "Explosion (then 
afterwards)". The top right panel shows that the GW amplitudes monotonically increase or decease with time in the 
"Explosion" phase reflecting the geometry of the expanding shock (bottom panel), which is either prolate (increase), 
oblate (decrease), or spherical (in-between). In each panel, the chosen progenitor models of It^I are given in each panel 
by 12, 20 Mq and the input luminosity is indicated by the value of Ly^ in unit of 10^^ erg/s. 
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shows matter GWs obtained for a 15 Mq progenitor model 1721] with neutrino luminosity of 



Ly^ = 3.7 X 10''^erg/s. Three typical GW emission cites are clearly presented, which are "Prompt 
convection", "Postshock convection" (unless it is shown in the plot), "Non-linear SASI with 
SASI plumes", and "Explosion". 

Shortly after the bounce shock is formed, negative entropy gradient behind the stalling shock 
predominantly gives rise to prompt convection. GWs indicated by "Prompt convection" (top left 
panel in Figure [T]) comes from this. Later on, the PNS convection driven by the negative lep- 
ton gradient near its surface and the neutrino-driven convection in the postshock heating region 
develop. This corresponds to the GW emission by "Postshock convection" (unless shown in 
the plot) after the prompt convection till ~ 550 ms after bounce. Subsequently, the GW ampli- 
tudes become much more larger as the SASI enters to a non-linear phase with violent sloshing 
of the postshock material. Large spikes appearing in the "Non-linear SASI" phase come from 
the down-flowing "SASI plumes" striking the PNS surface. Afterward when the sloshing shock 
turns to be explosion (indicated by "Explosion"), the sign of the GW amplitudes change reflect- 
ing the geometry of the expanding shock, that is either prolate (increase), oblate (decrease), or 
spherical (in-between) as shown in the top right and bottom panels in Figure [1] On top of these 
illuminating findings, they argued that the characteristic GW frequency has a tight correlation 
with a deceleration timescale of the accreting material at the lower boundary of the postshock 
convective region. As will be menti oned in section [J. L4I this trend is also observed in a recent 



3D simulation by Miiller et al. 113111 



3.1.2. Features of Neutrino GWs in experimental 2D simulations 

To understand the behavior of neutrino GW^ in 2D simulations, it is useful to look the stress 
formula. In 2D simulations, the nonvanishing component of the neutrino GWs can be expressed 



as um 



TT 4-G r , r , , diy{e,t) 



where G is the gravitational constant, c is the speed of light, R is the distance of the source to the 
observer, dly/dQ. represents the direction-dependent neutrino luminosity emitted per unit of solid 
angle into direction of £1, and 0(0 ) denotes the quantity, which depends on the angle measured 
from the symmetry axis (O ), 

0(6'') = ;rsin0'(-l +2|cos0'|). (2) 

It is easy to show that the neutrino GW is zero if the neutrino emission is isotropic [dly/ (0 , f )=constant). 
This function has positive values in the north polar cap for < ff < 60° and in the south polar 
cap for 120° < 0' < 180°, but becomes negative values between 60° < 0' < 120°. Therefore, 
the polar and equatorial excess in the neutrino emission tends to make a positive and negative 
change in the neutrino GWs, respectively. 

Figure |2] shows p roper ties of neutrino GWs obtained in the parametric 2D simulations by 
Kotake et al. (2009) 112611 who used the light-bulb scheme to obtain explosions. In their sim- 
ulations, the initial conditions were derived from a steady-state approximation of the postshock 



As anisotropic matter motions gen erate GWs , anisotropic neutrino e missio n also gives rise to GWs, which has been 
oiiginally pointed out in l ate 197Q's by Il32lll33ll (se e rece nt progress in Il34ll ). It is expected as a piimary GW source 
also in gamma-ray bursts QHIIli] and Pop m stars fBTll . 
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Figure 2: Gravitational waveforms obtained in 2D parametric explosion models in Kotake et at. (2009) 112611 (for a 
distance of 10 kpc). For the four models (of A, B. C, and D) with different input luminosities (see text for more details), 
the top panels show the GWs only from neutrinos (top left) and from the sum of neutrinos and matter (top right). The 
time is measured from the epoch when the simulation starts from the initial condition (see text). In all the models, the 
SASI enters to the non-linear regime at about 100 ms, simultaneously making the amplitudes deviate from zero. The 
bottom left panel shows the directional dependent neutrino luminosity: dly/dSl for model A, in which the polar excess 
from equator in the luminosity is indicated by "(South pole) - (Equator)" or "(North pole) - (Equator)". Vertical lines 
represent the epochs of ^ = 294 and 370 ms, respectively (see text for details). For the chosen two epochs, the bottom right 
panel shows temperature- (the left-half of each panel) and density- (the right-half) distributions in the meridian section. 
The insert of each panel shows the neutrino GW amplitudes, in which the green point indicates the time of the snapshot. 
The central region colored by white (50 km in radius) represents the inner boundary taken in the 2D simulation. 
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structure (not from stellar evolution calculations) and the dynamics only outside an inner bound- 
ary at 50 km was solved. To estimate the neutrino anisotropy (e.g., dly/{0 ,t) in Equation[T]i, a 
ray-tracing analysis was performed. In the two top panels of Figure |2] the input electron neu- 
trino luminosity {Ly} is indicated for four models A, B, C, and D, respectively corresponding to 
Lv„ = 6.4, 6.6, 6.7, and 6.8 x lO^^erg/s. 

The top left panel shows the GW amplitudes contributed only from neutrinos. Comparing the 
right panel, which shows the total amplitudes (neutrino + matter), it can be seen that the overall 
structures of the waveforms are predominantly determined by the neutrino GWs with slower 
temporal variations (> 50 ms), to which the GWs from matter motions with rapid temporal 
variations (< 10 ms) are superimposed. 

The bottom left panel shows the angle dependent neutrino luminosity (e.g., dly/dO. in Equa- 
tion ([T]i) for model A, in which the polar excess from equator in the luminosity is indicated by 
"(South pole) - (Equator)" or "(North pole) - (Equator)". As shown, the dominance of the neu- 
trino emission in the north (green line) and south poles (red line) are closely anti-correlated. This 
is the consequence of the low-mode nature of the SASI, here of ^ = 1 (in terms of an expansion 
in spherical harmonics with order {). In fact, the left-hand side in the bottom right panel shows 
that at 294 ms, the blob encompassing the regions inside the stalled shock is moving from the 
southern to the northern hemisphere, leading to the compression of the matter in the south hemi- 
sphere, which is vice versa at 370 ms (right-hand side). Recalling again that the function (Figure 
|2l) is positive near poles, the polar excess in the neutrino GWs makes the positively growing fea- 
ture depicted in the neutrino GWs (e.g., top left panel in Figure|2|. In addition, the waveforms 
show large negative growth for some epochs during the growth of SASI (e.g., for models C and 
D in the top left panel). Such a feature comes from a dominance of the neutrino emission in 
the equatorial direction^. Reflecting the nature of the SASI and convection both of which grow 
chaotically with time, they pointed out that there may be no systematic dependence of the input 
neutrino luminosities on the maximum GW amplitudes. 



3.1.3. GW signatures in 2D radiation-hydwdynamic simulations 

Having summarized the findings in parametric 2D simulations, we are now ready to discuss 
GW signatures obtained in 2D radiation-hydrodynamic simulations in which a spectral neutrino 
transport is solved. Similar to Table 1, but Table 2 shows a summary for some selected GW 
predictions obtained in the first-principle 2D model0 Since the model predictions are not too 
many (as of 201 1), the important GW features obtained mostly by the MPA simulation and also 
by the OakRidge+ simulation are all summarized in this section by reproducing their results as 
much as possible. 

The first example is for a rotating I5M0 progenitor model by the MPA simulation (e.g., 
¥rh . see the numerical details in Table 1), the detailed GW analysis of which was performed 
by Miiller et al. (2004) (Figure O. The initial angular velocity for the model was taken to be 
0.5 rad/s, which is slightly faster (a factor 5) than the one predicted in a recent stellar evolution 
calculation fsJl- As shown in the top left panel, they found a quasi-monotonically increasing 
GW signal at the end of their simulation (~ 270 ms after bounce), which predominantly comes 
from neutrino GWs. I n the case of r otational core-collapse, neutrino emission can be stronger 
along polar directions 1 138 , 139l 140ll . As discussed in section [J.I.2l the polar excess in neutrino 



To capture this trend, the ray-tracing calculation would be more accurate than a ray-by-ray approach. 
"Note that much more stronger GW emission than those in Table 2 was expected if the acoustic mechanism would 
work. See section 7.1 in j22ll for more details. 
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Progenitor 


Group 


Mechanism 


Dim. 


ffin 


|/2maxl(cm) 


/peak 


(rotation) 


(Year) 


(Domain) 


(Hydro) 


(ms) 


(from) 


(Hz) 


11.2Mo^ 


MPA[771 


v-driven 


2D 


-200 


-30 


-700 


(No) 


(2004) 


(180°wedge) 


(PN) 




(matter) 






Oak Ridge+ 


v-driven 


2D 


-550 


80 


-900 


(No) 


[801(2009) 


(180°wedge) 


(PN) 




(neutrino) 




13 


NAOJ+ 


v-driven 


2D 


-700 


18 


-600 


(Both) 


[8]J(2011) 


(180° wedge) 


(N) 




(neutrino) 




15Mo|73J 


MPA[77] 


v-driven 


2D 


-270 


-150 


-600 


(Yes) 


(2004) 


(90°wedge) 


(PN) 




(neutrino) 






MPAITTJ 


v-driven 


2D 


-400 


-130 


-600 


(No) 


(2006) 


(180°wedge) 


(PN) 




(neutrino) 




15 MoQ2] 


Oak Ridge+ 


v-driven 


2D 


-550 


>200 


-900 


(No) 


[801(2009) 


(180° wedge) 


(PN) 




(neutrino) 






Oak Ridge+ 


v-driven 


2D 


-550 


-200 


-900 


(No) 


[80] (2009) 


(180° wedge) 


(PN) 




(matter) 





Table 2: Similar to Table 1. but for selected GW predictions obtained in the state-of-the-art 2D neutrino-radiation hy- 
drodynamic simulations. See Table 1 for the details of the employed numerical techniques. The column of (Domain) 
indicates that the computational volume in the lateral direction is solved fully < 9 < 180° ("180°wedge") or with an 
assumption of certain imposed symmetry ("90°wedge"). ffi„ represents the final simulation time, and |/imaxl indicates 
the absolute maximum GW amplitudes which come either from "neutrino" or "matter" GWs. /peak represents the peak 
frequency in the GW spectra. 



14 



emission could be a possible explanation of the increasing trend. It is also worth mentioning that 
this f eature can be seen in the waveform for a rapidly-rotating accretion-induced collapse model 
1 14111 of the "Princeton-H" simulation (e.g., Figure 10 in 12211 ). 

The top right panel of Figure [3] shows the growth of large scale, non-radial pulsations in the 
post-shock region (bright regions in the entropy plot (600x600 km)), which leads to a slower 
temporal evolution (> 10 ms) of the neutrino GWs (see the line labeled by "v's" in the top left 
panel). As seen, an equatorial symmetry was imposed for the 2D simulatior^, which is indi- 
cated by "90°wedge" in Table 2. Comparing the GW spectra from the sum of matter motions 
and neutrinos (bottom left panel in Figure O to the one only from neutrinos (bottom right), one 
can clearly see that the low-frequency part of the spectrum (below ~ 100 Hz) is dominated by 
the contribution of the neutrino GWs. As mentioned above, this is because their wave amplitude 
varies on much longer timescales (> 10 ms) than that of the mass flow (~ ms, e.g., the wave- 
form in the top left panel). For this model, the maximum amplitude comes from the neutrino 
contribution (e.g., the top left panel and |/imaxl with the comment of (neutrino) in Table 2). 

It should be noted whether the maximum amplitude comes from matter or neutrino contri- 
bution change from models to models as shown in Table 2 (see the sixth column indicated by 
(from)), which may reflect the stochastic nature of the SASI and convection which is determined 
by the non-linear hydrodynamics. In contrast, the peak frequencies in the GW spectrum (typ- 
ically in the range of 600-1000 Hz) always come from the matter GWs (albeit not explicitly 
mentioned in the table) due to its short temporal variations (/peak in Table 2). The bottom panels 
of Figure [3] show that while the signal-to-noise ratio is probably too small for this event to be 
detectable by LIGO I, it could be well detected by the Advanced LIGO for a Galactic source. 
This is a generic feature of the detectability of the GWs from neutrino-driven explosion models. 

Figure |4] is the same as Figure [3] but for a non-rotating exploding model of an 11 .2 Mq 
progenitor of Buras et al.|l77|]. In contrast to the rotating model mentioned above, no signal is 
emitted near at bounce for this non-rotating model (compare the insert in the top left of Figure[3] 
and the waveform near f = in the top left panel of Figure|4|. For this light progenitor model, the 
wave amplitude from prompt convection (till ~ 50 ms after bounce) is very small. Only after the 
neutrino-driven convection becomes stronger (~ 80 ms after bounce), the maximum (absolute) 
amplitudes rise to the values of several 10 cm. Although the slower temporal evolution of the 
neutrino GWs is similar to the one in Figure [3j the monotonically growing trend is not seen for 
this model. The GW spectrum (bottom panels in Figure HJi is qualitatively quite similar to that 
of the rotating model (FigureO, but the GW amplitudes become much smaller. This is probably 
because of a smaller core mass of the 11.2 M© progenitor, which could potentially make the mass 
quadrupole and the emergent neutrino luminosity smaller (compare the spectrum of neutrino GW 
in the bottom right panel of Figure [3] and |4|. 

The left and middle panels in Figure |5] show the wave amplitude (blue line (total), red line 
(only from neutrinos)) for a non-rotating exploding 15 Mq progenitor in Marek & Janka (2006) 
18211 in which a soft variant of the Lattimer-Swesty EOS (left) and a rather stiff Hillebrant- Wolff 
EOS (middle) was employed|£3 Although an overall trend is similar with each other, they found 
several important EOS effects on the GW signatures. The first one is about the GWs from prompt 
convection (e.g., the signals before 50 ms after bounce). As seen, the wave amplitude for the 
stiffer EOS (~ 20 cm) becomes up to a factor of 2 larger maximum amplitudes than its softer EOS 



^"Hence the dynamics is nortli-south symmetric. 

^ ' An incompressibility at nuclear densities (A") is 1 80 MeV and 281 MeV for the Lattimer-Swesty and Hillebrant- Wolflf 
EOS, respectively. 
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Figure 3: Gravitational waveforms from neutrinos (indicated by "v's" in tlie top left panel) or from matter motions 
(indicated by"flow" in the panel) and their spectra (sum of matter and neutrino (bottom left) and contri butio n only from 
neutrino GWs (bottom right)) for a rotating model (sl5r) of 15 Mq progenitor in Miiller et al. (2004) 112011 . The insert 
in the top left panel shows an enlargement of the signal near bounce, which comes from prompt convection (e.g., section 
13.1. It . Top right panel shows the growth of large scale, non-radial pulsations in the post-shock region (bright regions in 
the entropy plot (600x600 km)), which leads to a s lowe r temporal evolution (> 10 ms) of the neutrino GWs. All of these 
plots are by courtesy of Miiller and the coauthors Il20ll . 



counterpart (~ 10 cm). They pointed out that the prompt convection exhibits larger amplitude, 
slightly higher frequencies, and more power for the stiffer EOS. They found that this is because 
the region of prompt convection involves more mass and extends to larger radii. This may agree 
with one's intuition that the inner-core mass near bounce becomes larger in the case of stiffer 
EOS (due to the suppressed electron capture), potentially providing a larger quadrupole GW 
emission. 

The maximum (absolute) amplitudes come from the neutrino contribution irrespective of 
the employed EOSs (see the column of (from) for the non-rotating 15Mq progenitor in Table 
2). The maximum (absolute) amplitude is about 130 and 60 cm for the softer and stiffer EOS, 
respectivelyl^ The higher neutrino luminosities that are radiated from a more compact PNS 
in the softer EOS predominantly give a favourable condition for a more efficient emission of 
the neutrino GWs. However, their peak frequencies are typically below ~ 100 Hz (e.g., the 
bottom right panel in Figure |3]l, making the component of the neutrino GWs very difficult to be 
detectable for ground-based detectors whose sensitivity is limited mainly by seismic noises at low 



■Note that the result only for the softer EOS is given in Table 2. 
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Figure 4: Similar to Figure [3] bu t for a non-rotating explosion model in the MPA simulations of an 11.2 M© progenitor 
(taken from Miiller et al. (2004)fT20l|). The top right panel indicates the success explosion of this model, in which the 
entropy (bright) with different timescales (181 (top) or 225 ms (bottom) after b ounc e) are displayed. The polar axis is 
directed horizontally. These figures are by courtesy of Miiller and the coauthors I120II . 



frequencies. Therefore the matter GW signals are more important in discussing the detectability. 
For a Galactic source, the matter signal are of the orders of 10"^^ regardless of the different 
EOSs. They discussed that they should be marginally observable with LIGO 1 but it is not easy 
to tell the EOS difference. This is because although the spectral peak become (slightly) higher if 
the EOS is soft, this simultaneously sticks out of the highest sensitivity frequency domain of the 
LIGO instrument. 

In contrast to an anticipation in section U. 1.21 their neutrino GWs decrease with time. They 
argued that this comes from the excess of neutrino emission in the equatorial direction, which 
was argued to be determined by the anisotropic transport of heavy-lepton neutrinos in the vicinity 
of the PNS E§. To pin down the GW signal from neutrinos, one apparently needs to precisely de- 
termine the neutrino anisotropy which should require multi-angle neutrino transport calculation. 
Currently it would be practicable in 2D simulations Isill if enough number of momentum-space 
angles could be cast (probably) by using the next-generation supercomputers. But before that, 
we may need to be a little bit careful in interpreting the neutrino GWs. 



On the other hand, they did see the stronger emission in the polar regions regarding v,, and v,,, which is consistent 
with the discussion in sec tion [3. 1.21 
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Figure 5: Gravitational wavefomis for a non-rotating exploding model in the MPA simulations of a 15 Mq progenitor 
1831 in whic h the Lattimer-Swesty (left panel) or Hillebrant-WolfF EOS (middle panel) was employed (taken from Marek 
et al. (2009) Il28ll ). The right panel indicates the explosion of this model, in which the lefthand panel shows color-coded 
entropy distribution, the righthand panel the radial velocity components with white and whitish hues denoting matter at 
or near rest, w hile black an'ows indicate the direction of the velocity fields. These figures are by courtesy of Mai'ek and 
the coauthors Il28ll . 



Keeping this in mind, the neutrino GWs (blue line in the left panel of Figure |6]l for a non- 
rotating 15Mq exploding model in the OakRidge-i- simulation (e.g.. Table 1 for the numerical 
details) show again a monotonically increasing trend. They obtained a similar trend in the wave- 
forms of their exploding models of 12 and 25 Mq progenitors |f89'l. In consistent with the analysis 
in section lT. 1.11 the matter GWs show a monotonic increase after the S ASI-driven sloshing shock 
turns into 2D prolate explosions (> 300 ms postbounce in the left panel of Figure |6|l. For the 
first time this model succeeded in capturing the feature, which was not seen in other radiation- 
hydrodynamic simulations. To produce the increasing trend, the anisotropic kinetic energy of the 
expanding shock needs to be enough large. Relatively earlier shock revival with larger explo- 
sion energies (e.g.. Table 1) obtained in the OakRidge-i- simulations could satisfy the condition. 
The maximum amplitudes for a Galactic supernova reach to the orders of 10"^' regardless of the 
employed progenitors, which could be a promising target for the advanced LIGO (bottom right 
panel in Figure |6]l. 

Figure [7] shows the wave amplitudes (left panels) with their spectra (middle panels) for a 
non-rotating (top panels) and rapidly rotating model (bottom panels) for a 13 Mq progenitor 117 111 
obtained in the NAOJ-i- simulation (e.g.. Table 1 for the numerical details). For the non-rotating 
model, the neutrino GWs show the increasing trend (green line in the top left panel). This model 
explodes in a unipolar manner (top right panel), expelling material only in the northern part of 
the core. Due to this one-sided explosion with smaller mass ejection, this model explodes only 
weakly with decreasing explosion energies less than ~ 4 x 10'*'' erg at the final simulation time 
(see Figure 4 in jsiH ). This should be the reason why the matter GW signals (red Une in the top 
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Figure 6: Similar to Figure|3]but for a non-rotating exploding model in the OakRidge+ simulations of a non-rotating 15 
Mq progenitor (from Yakunin et al. (2009) (s^ ). The left panel shows the gravitational waveforms contributed from 
matter motions (indicated by matter) and from neutrinos (indicated by neutrino). The insert shows an enlargement of 
the signal near bounce, which comes from prompt-convection. The top right panel shows entropy distribution at 244 ms 
after bounce for the progenitor. A bipolar explosion with high-entropy outflows (yellow-orange-red) is shown along the 
symmetry axis. The bottom right panel shows the spectra with the sensitivity curve of the advanced LIGO. These figures 
are by courtesy of Yakunin and the coauthors i89h . 



left panel) are much smaller than those obtained in the Oak Ridge-H simulation (e.g., red line in 
Figure |6]l. 

On the other hand, a bipolar explosion is seen to be obtained for their 2D model with rapid 
rotatioij^ (bottom right panel in Figure |7]i. Since the north-south symmetric {{ - 2) explosion 
can expel more material than for the unipolar explosion, the explosion energies becomes larger 
(e.g.. Figure 4 in ISlll ). simultaneously making the matter GWs larger In fact, the increasing 
trend in the matter GWs is seen after ~ 400 ms after bounce (red line in the bottom left panel in 
Figure|7]i. The neutrino GWs firstly decrease after bounce (blue line in the same panel), but later 
shift to exhibit the increasing trend (> 500 ms after bounce). The middle panels (for a supernova 
at a distance of 10 kpc) show that the spectrum rises to a broad peak between ~ 600 and 900 
Hz, which cou ld be visible to the advanced-class detectors for the galactic source (see further 
discussions in il42ll ). 



The precollapse angular velocity was taken to be 2 rad/s. 
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Figure 7: Gravitational waveforms (left panels) with their spectra (middle panels) for a non-rotating (top panels) and 
rapidly rotating model (bottom panels) for a 13 Mq progenitor ItTIi obtained in the NAOJ+ simulation (e.g., Table 1). 
For the rotating model, the precollapse angular velocity is taken to be 2 rad/s. Right panels show snapshots of the density 
(left half) and the entropy (right half) for models without (top) or with rotation (bottom) at the epoch when the shock 
reaches 1000 km, corresponding to ~ 470 ms after bounce in both cases. The m iddle panels are for a supernova at a 
distance of 10 kpc. These figures are by courtesy of Suwa and the coauthors Il42ll . 



3.1.4. Stochastic GW nature in experiment 3D simulations 

(By definition) 2D simulations so far mentioned in previous sections have an assumption of 
axisymmetry. Due to the symmetry axis, the growth of the SASI and the large-scale postshock 
convection could develop along the axis preferentially, thus suppressing anisotropics in explo- 
sions. Since the GW emission is very sensitive to the degree of the explosion anisotropics, 3D 
simulations are apparently needed for a more accurate GW prediction. 

Since it is still computationally very expensive to solve the spectral neutrino transport, a 
Ught-bulb scheme (e.g., ll54ll ) has been often employed in 3D simulations to trigger explosions, 
in which the heating and cooli ng by neutrinos are treated by a parametric manner (e.g., section[2] 
and references in liiliillil). Using the recipe, a handful of GW predictions based on the 3D 
simulations have been reported so far, which we are going to review in this section. As a preface, 
we just hke to mention that although these studies use experimental approaches and thus their 
findings should be tested by the first-principle 3D simulations in the future (e.g., ll79ll ). they have 
indeed shed a new light towards a better understanding of the GW signatures as they have done 
so in 2D simulations (see disc ussio ns in sections D . 1 . II and |3 . 1.2b . 

Miiller and Janka (1997) 111 1711 coined the first study to analyze the GW signature of 3D 
non-radial matter motion and anisotropic neutrino emission from prompt convection in the outer 
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layers of a PNS during the first 30 ms after bounce. Their first 3D calculations using the light- 
bulb recipe were forced to be performed in a wedge of opening angle of 60°. Albeit with this 
limitation (probably coming from the computer power at that time), they obtained important 
findings that because of smaller convective activities inside the cone with slower overturn veloc- 
ities, the GW amplitudes of their 3D models are more than a factor of 10 smaller than those of 
the corresponding 2D models, and the wave amplitudes from neutrinos are a factor of 10 larger 
than those due to no n-radial matter motions. With another pioneering (2D) study by Burrows & 
Hayes (1996) jl 1811 . it is worth mentioning that those early studies had brought new blood into 
the conventional GW predictions, which illuminated the importance of the theoretical prediction 
of the neutrino GWs. 



A series of findings obtained by Fryer et al. in early 2000s 11551114411 have illuminated also 



the importance of the 3D modeling. By running their 3D Newtonian Smoothed Particle Hy- 
drodynamic (SPH) code coupled to a gray flux-limited neutrino transport scheme, they studied 
the GW emission due to the inhomogeneous core-collapse, core rotation, low-modes convection, 
and anisotropic neutrino emission. Although the early shock-revival and the subsequent pow- 
erful explosions obtained in these SPH simulations have yet to be confirmed by other groups, 
their approach paying particular attention to the multiple interplay between the explosion dy- 
namics, the GW signatures, the kick and spins of pulsars, and also the non-spherical explosive 
nucleosynthesis, blazed a new path on which current supernova studies are progressing. 

We also st udied the GW s ignals from 3D models that mimic neutrino-driven explosions aided 
by the SASI (iM 11271 114511 . These studies were an extension of the 2D simulations already 
mentioned in section 13.1.21 in which the light-bulb scheme was used to obtain explosions and 
the initial conditions were derived from a steady-state approximation of the postshock structure 
and the dynamics only outside an inner boundary at 50 km was solved. The left panel of Figure 
[8] shows the evolution of 3D hydrodynamic features from the onset of the non-linear regime of 
SASI (top left) until the shock break-out (bottom right) with the gravitational waveform from 
neutrinos inserted in each panel. As seen, the major axis of the growth of SASI is shown to be 
not aligned with the symmetric axis (:Z axis in the figure) and the flow inside the standing shock 
wave is not symmetric with respect to this major axis (see the first and third quadrant). This 
is a generic feature obtained in their 3D models, which is in contrast to the axisymmetric case. 
As discussed in section [3. 1 .21 the GW amplitudes from the SASI in 2D showed an increasing 
trend with time due to the symmetry axis, along which SASI can develop preferentially. Free 
from such a restriction, a variety of the waveforms is shown to appear (see waveforms inserted 
in Figure[8ll. Furthermore, the 3D standing shock can also oscillate in all directions, which leads 
to the smaller explosion anisotropy than 2D. With these two factors, the maximum amplitudes 
seen either from the equator or the pole becomes smaller than 2D. On the other hand, their sum in 
terms of the total radiated energy are found to be almost comparable between 2D and 3D models, 
which is likely to imply the energy equipartition with respect to the spatial dimensions. 

The right panel of Figure [8] shows the gravitational waveforms for models with different 
neutrino luminosity. The input luminosity for the bottom panel is smaller than that for the top 
panel only by 0.5%. Despite the slight difference, the waveforms of each polarization are shown 
to exhibit no systematic similarity when seen either from the pole or equator. 

The stochastic nature of GWs produced by the interplay of the SASI and neutrino-driven 
explosions have been also confirmed more recently by Miiller et al. (201 1) il3lil . They ana lyzed 
the GW signatures based on their 3D models of parametric neutrino-driven explosions 11 14611 
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Figure 8: Left panel shows four snapshots of the e ntrop y distributions of a representative 3D supernova explosion model 
(con'esponding to model A in Kotake et al. (2009) |127tl ). The second and fourth quadrant of each panel shows the surface 
of the standing shock wave. In the first and third quadrant, the profiles of the high entropy bubbles (colored by red) inside 
the section cut by the ZX plane are shown. The side length of each plot is 1000km. The insets show the gravitational 
waveforms from anisotropic neutrino emissions, with '+' on each curves representing the time of the snapshot. Note that 
the colors of the curves ai'e taken to be the same as the top panel. After about 100 ms, the deformation of the standing 
shock becomes remarkable marking the epoch when the SASI enters the non-linear regime. At the same time, the 
gravitational amplitudes begin to deviate from zero. Right panel shows gravitational waveforms from neutrinos (bottom) 
and from the sum of neutrinos and matter motions (top), seen from the polar axis and along the equator (ind icated by 
'Pole' and 'Equator') with polarization (+ or X modes) for two representative 3D models of A and B (see Il27ll for 
details), in which the input luminosity for the two pair panels differs only 0.5%. From the right panels, it can be seen 
that the overall structures of the waveforms are predominantly determined by the neutrino GWs which was also the case 
in 2D simulations (e.g., section [3. 1.21 . The distance to the source is assumed to be 10 kpc. 
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in which the gray spectral neutrino transport is solved in a ray-by-ray manner 15911^^1 . At the 
sacrifice of cutting out the high-density core (at around 80 - 120 km in radius), their 3D models 
succeeded in following the dynamics from the post-bounce accretion phase through the onset of 
the explosion more than one second for the first time. Among a number of interesting findings 
obtained in their work, we here focu s on the progenitor dependence that could not be studied by 
Kotake et al. (2009,201 1) il27lll45ll who used the ideahzed initial conditions. 

The top two panels in Figure |9] show snapshots of blast morphologies between two models, 
in which a non-rotating 15 Mq progenitors either evolved by Woosley and Weaver (1995) (left) 
or by Limongi et al. (1995) (right) was employed, respectively. Note that they succeeded in 
following an unprecedentedly long-term 3D evolution (up to 1 .4 s aft er bounce) by utilizing 
an axis-free coordinate system; the so-called Yin- Yang grids (see lll46ll for more details). The 
middle panels in Figure |9] indicate that the total GW amplitude is significantly different from 
that of the flow-only GW amplitudes (not shown) for models which do or do not exhibit PNS 
convection below the neutrino sphere. Particularly at late times (f > 0.8 s in the middle panels), 
anisotropic neutrino emission causes a continuing growth of the GW amplitudes for the model 
using the Woosley-Weaver (indicated by W15) progenitor (left) instead of saturation seen for the 
model using the Limongi progenitor ( indicated by L15, right panel). By a detailed spectrogram 
analysis of the GW energy distributions (see their Figure 12), they revealed that the primary agent 
to make the model discrepancy is the PNS convection, which develops much more vigorously for 
the W15 progenitor model than for the L15 progenitor model. This can be also depicted in the 
bottom panels of Figure |9] which shows bigger neutrino anisotropy for the W 1 5 progenitor model 
(left) in the late postbounce phase (f > 0.8 s) than for the L15 progenitor model (right). They 
furthermore pointed out that in their 3D models, very prominent, quasi-periodic sloshing motions 
due to the SASI are absent and the emission from different surface areas facing an observer adds 
up incoherently, so that the measurable modulation amplitudes of the GW signals (as well as 
neutrinos) are significantly smaller than predicted by 2D simulations ( I129[l89iil ). 

The effects of stellar rotation on the stochastic nature of the GWs have been recently studied 
by Kotake et al. (2011) il45ll who used the same numerical techniques that was already men- 
tioned in explaining Figure |8] The top two panels in Figure [TO] show the gravitational waveform 
for a typical 3D with rotation (left:total amplitudes, right:neutrino only). To construct a model 
with rotation, a uniform rotation was manually given on the flow advecting from the outer bound- 
ary of the iron core as in 11671] . whose specific angular momentum is assumed to agree with recent 
stellar evolution models ll84ll . Comparing to the right panel in Figure [8] (i.e., in the absence of 
rotation), one can clearly see a sudden rise in the GW amplitude after around 500 ms for the 
rotating model (blue line in the top panels in Figure [TOb . which is plus mode of the neutrino 
GWs seen from the equator. These features were found to be common in the fifteen 3D models 
computed in their study. 

Seen from the equatorial direction (the middle left panel of Figure [TOl i. one may guess the 
presence of the sloshing modes, but it just happens to develop along the rotational axis (z-axis) 
at this epoch. It should be emphasized that the dominance of /i^''" observed in the current 3D 
simulations have nothing to do with the one discussed in section 13.1.21 Free from the 2D axis 
effects, the major axis of the SASI changes stochastically with time, and the flow patters behind 
the standing shock also change in every direction. As a result, the sloshing modes can make only 
a small contribution to the GW emission. The remaining possibility is that the spiral flows seen 



^^It is worth mentioning tliat a complete derivation of formulae for extracting th e neu trino GWs is given for arbitrary 
positioning between the source and observer coordinate systems (e.g., section 4 in Il3lll ^. 
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Figure 9: Top two panels show snapshots illustrating a globally asymmetric 3D explosion in the post-accretion phase of 
two different progenitor models u sing the Woosley- Weaver (left) or Limongi progenitor (right) of a 15 Mq prog enitor 
(taken from Miiller et al. (2011) Il3lll ). Note that "-4" or "-3" differ only by the initial seed perturbations (see Il3lll 
for more details). The time and linear scale are indicated in each plot. Middle panels show the plus (blue) or cross (red) 
mode of the wave amplitude due to anisotropic mass flow and neutrino emission for models W15-4 (left) and L15-3 
(right), respectively. The solid curves show the amplitudes for an observer above the north pole of the source, while the 
other curves give the amplitudes at the equator. Similar to the middle panels, the bottom panels shows the anisotropy 
param eter of neutrino emission (only the plus mode is shown). These figures are by courtesy of Miiller and the coauthors 
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Figure 10: Top panels show gravitational waveforms from the sum of neutrinos and matter motions (left) and onl y fro m 
neutrinos (right) obtained in a parametric explosion 3D model with rotation (taken from Kotake et al. (2011) Il45ll ). 
Similar to the left panel in Figure[8] the middle panels show a partial cutaway of the entropy isosurfaces and the velocity 
vectors on the cutting plane for the equatorial (left) and polar observer (light), respectively. The bottom panel shows 
spectral distributions for t he m odel from matter motions ("Matter") and neutrino emission ("Neutrino") seen from the 
pole or the equator (e.g., Il45ll ) with the expected detec tion limits of TAMA300 flTll . first LIGO and advanced LIGO 
|147|| . LCGT ffsll . and Fabry-Perot type DECIGO Il48ll . The distance to the supernova is assumed to be 10 kpc. Note 
that for the matter signal, the + mode seen from the polar direction is plotted. 
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in the middle right panel should be a key importance to understand the GW feature mentioned 
above. In fact, by analyzing the matter distribution on the equatorial plane, it was found that the 
compression of matter is more enhanced in the vicinity of the equatorial plane due to the growth 
of the spiral SASI modes, leading to the formation of the spiral flows circulating around the spin 
axis with higher temperatures. As a result, the neutrino emission seen parallel to the spin axis 
becomes higher than the ones seen from the other direction. Remembering that the lateral-angle 
(ff) dependent function of the GW formulae (e.g., in Equation ([T]i) is positive near the north and 
south polar caps, the dominance of the polar neutrino luminosities leads to make the positively 
growing feature of /i^''" in the top panels of Figure[TO](blue line). From the spectral analysis of the 
gravitational waveform (the bottom panel of Figure [TOb. it can be readily seen that it is not easy 
to detect these neutrino-originated GW signatures with slower temporal evolution (> (9(10)ms) 
by ground-based detectors whose sensitivity is limited mainly by the seismic noises at such 
lower frequencies. However these signals may be dete ctabl e by the recently proposed future 



space interferometers like Fabry-Perot type DECIGO ( 1114811 . black line in the bottom panel). 
Contributed by the neutrino GWs in the lower frequency domains, the total GW spectrum tends 
to become rather flat over a broad frequency range below ~ 100 Hz. These GW features obtained 
in the context of the S ASI-aided neutrino-driven mechanism are different from the ones expected 
in the other candidate supern ova rn echanism, such as the MHD mechanism (e.g., section [X2] | 



and the acoustic mechanism ( 1112311 , e.g., section 7.1 in 112211 ). Therefore the detection of such 
signals could be expected to provide an important probe into the explosion mechanism. 

Finally, most of the 3D models summarized in this section cut out the PNS and the neutrino 



trans port is approximated by a simple light-bulb scheme 11 14511 or by the gray transport scheme 



1 13111 . Needless to say, these exploratory approaches are but the very first step to model the 
neutrino-heating explosion and to study the resulting GWs. As already mentioned, the excision 
of the central regions inside PNSs truncates the feedback between the mass accretion to the PNS 
and the resulting neutrino luminosity, which should affect the featu res o f the neutrino GWs. By 



the cut-out, eflicient GW e mission of the oscillating neutron star 1112311 and non-axisymmetric 



instabihties il22lllllllll2ll of the PNSs, and the enhanced neutrino emissions inside the PNSs 



1 12811 cannot be treated in principle. To elucidate the GW signatures in a much more quantitative 
manner, full 3D simulations with a spectral neutrino transport are apparently needed. This is 
unquestionably a vast virgin territory awaited to be explored for the future. 

3.2. GWs from MHD Explosions 

As already mentioned in the beginning of section |3] state-of-the-art results concerning the 
GW signatures emitted near bounce in the case of rapidly rotating core-collapse have been given 
in the recent review by Ott (2009) |5Ml (e.g., section 4). Conventionally the waveforms of the 
bounce signals are categorized into the three types, namely types I, II, and III. As explained in 
Ott (2009) H^j, type II and III waveforms are shown less likely to appear than type I, because 
a combination of general relativi ty and electron capture near core bounce suppresses multiple 



bounce in the type II waveforms 10711 . In general, a realistic nuclear equation of state 



(EOS) is stiff enough to forbid the type III waveforms. So the generic type of the bounce signals 
is now known to take the type I waveform. 

For avoiding overlap, we mainly focus on the postbounce GW emission in the context of 
rapidly rotational core-collapse in this section. The major GW emission cites have been proposed 
to be magnetohydrodynamic (MHD) outflows and nonaxisymmetric instabilities, which will be 
separately described in sections D . 2 . II and 13.2.21 respectively. 



26 



3.2.1. MHD outflows 

Numerical simulations of MHD stellar ex plos ions have a long histor\P^ started alrea dy in 
early 197 0's by LeBlanc and Wilson (1970) 114911 shortly after the discovery of pulsars il50l 



151L 115211 . However, it is rather only recently that the MHD studies come back to the front- 
end topics in the supernova research followed by a number of extensive MHD sim ulations (e.g., 



jp erno ' 

liMIiilliMImliMIiiSSIiiSI^ for references 



therein). Main reasons for this activity are o bservation s indicating very asymmetric explosion s 
lEs . 261. and the interpretation of magnetars 1 167 . 168 1 and gamma-ray bursts (e.g., 1 1 14 , 1 15 11 
as a possible outcome of the magnetorotational core-collapse of massive stars. 

The MHD mechanism of stellar explosions relies on the extraction of rotational free energy 
of collapsing progenitor core via magnetic field s. H ence a high angular momentum of the core is 
preconditioned for facilitating the mechanism lll69[l . Given (a rapid) rotation of the precollapse 
core, there are at least two ways to amplify the initial magnetic fields to a dynamically important 
strength, namely by the field wrapping by means of differential rotation t hat n aturally develops 
in the collapsing core, and by the magnetorotational instability (MRI, see lll70ll ). 

In the case of canonical initial magnetic fields (~ 10^G) l84ll . the fastest grow ing modes of the 
MRI are estimated to be at most several meters in the collapsing iron core 117111 . At present, it is 
generally computationally too expensive to resolve those small scales in the global MHD simula- 
tions, typically more than two or three orders-of-magnitudes smaller than their typical finest grid 
size. Let us remind first that the MHD explosions presented in the rest of this section are predom- 
inantly generated by the field wrapping mechanism assuming a very strong precollapse magnetic 
field (Bo ^ 10" G). Recent stellar evolution calculations show that this e xtreme condition 
could be really the case, albeit minor (~ 1% of massive star population 117211 ). for progenitors 
of rapidly rotating m etal-poor stars, which experience the so-called chemically homogeneous 
evolution ifTTl [Tlal . 

Figures[TT|to[T3]show several examples of gravitational waveforms obtained in representative 
MHD explosion models, in which the left and right panel shows the waveform and the blast 
morphology, respectively. After the bounce and ring-down signals with their typically duration 
of ~ 10 ms after bounce, a quasi-monotonically growing trend can be seen in the every waveform 
in the left panels of Figure [TTItofTS] 

Figures [TT] and [12] are from Obergaulinger et al. (2006) and Shibata et al. (2006) who 
obtained 2D MHD explosions either in Newtonian and approximate GR or full GRMHD simula- 
tions, respectively. Both of them employed polytropic precollapse models, a phenomenological 
EOS that mimics deleptonization and neutrino cooling. They pointed out that the increasing trend 
comes from bipolar flows driven by MHD explosions, which can be visible only for cores with 
preco llapse magnetic fields over Bo ^ 10'^ G. 3D MHD simulations by Scheidegger et al. (2010) 



1 12 1 included realistic EOSs and a deleptonization effect based on ID-Boltzmann simulations 



1 1311 . Compared to the corresponding 2D models, they pointed out that the jet-like explosions 



in their 3D models are much more difficult to obtain because the wind-up of the poloidal into the 
toroidal field does not proceed efficiently enough due to the growth of spiral SASI modes. But 
when a very rapid precollapse angular velocity of fio - 37r rad/s and (strong) magnetic field of 
Bo = 10'~ G were assumed, they observed the secularly growing trend associated with the MHD 
explosions also in their 3D simulations (Figure [T3]). 

In most of the mentioned MHD simulations, the bounce shock generally does not stall and 
a prompt explosion occurs within a few ten milliseconds after bounce. In fact, the typical 



'almost as old as the neutrino-driven model in Colagate and White (1966). 
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Figure 11: L eft panel shows a gravitational waveform in the 2D post-Newtonian simulations by Obergaulinger et al. 
(2006) I159II for their model A1B3G3-D3M13 in terms of the quadrupole amplitude Af^^. The dashed line comes from 
the contiibution of regions inside 60 km in radius while the soUd line shows the total signal. Comparing to the right panel 
(velocity contours at ~ 50 postbounce), the postbounce offset in the waveform is shown to come from the outer regions 
with bipolar outflows. These figures are by courtesy of Obergaulinger and the coauthors. 




Figure 12: Left panel shows a gravitational waveform in the 2D GRMHD simulations by Shibata et al. (2006) Il73ll for 
their models Al to A4. For models with stronger precollapse magnetic fields (pink and blue lines), the quasi-increasing 
trend can be clearly seen. Right panel shows density contours with velocity vectors for a snapshot (t = 152 ms for model 
A4) illustrating MHD-driven jet-like outflows. These figures are by courtesy of Shibata and the coauthors. 
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Figure 13: Left pair panels show gravitational waveform in the 3D post-Newtonian simulations by Scheidegger et al. 
2010) fTT2h for their model R4E1CF in which the left- and right- handside corresponds to the total amplitude and contri- 
bution only from the magnetic field, respectively. The right panel shows the 3D blast morphology of thek MHD explosion 
(entropy distribution at ~18 ms postbounce for a cubic of 300^ km^). These figures are by courtesy of Scheidegger and 
the coauthors. 



timescales shown in Figures [TTI to [T3] ar e at most 20 ms after bounce. This is because some 
of them followed adiabatic core-collapse lll59Ul73ll . in which a p olytro pic EOS is employed to 
mimic supernova microphysics, and the deleptonization scheme in 111 1211 cannot capture a drastic 
energy-loss behind the bounce shock due to neutrino cooling. However, for models with weaker 
precollapse magnetic fields akin to the current GRB progenitors, the prompt shocks stall firstly 
in the core like a conventional supernova model with more soph istica ted neutrino treatment (e.g., 
the MHD version of the Princeton-t- simulations presented in lll60ll '). In such a case, the onset 
of MH D explosio ns, depending on the initial rotation rates, can be delayed till ~ 100 ms after 

bounce ifieafiel . 



Takiwaki and Kotake (201 1) 017411 addressed this issue by performing 2D special relativistic 
MHD (SRMHD) simulations with the use of an approximate GR potential 1115811 in which a 
neutrino leakage scheme was employed to mimic neutrino cooling. The top left panel of Figure 
[14] shows the quasi-monotonically increasing trend, which is obtained for a model with strong 
precollapse magnetic field {Bq = lO'^G) also with rapid rotation initially imposed (J3 parameter = 
0.1% with /3 representing ratio of the rotational energy to the absolute value of the gravitational 
energy prior to core-collapse). Such a feature cannot be observed for a weakly magnetized model 
(Bo = 10"G. top right panel). 

To understand the origin of the increasing trend, it is most straightforward to look into the 
quadrupole GW formula, which can be expressed in 2D as 117511 . 



h - — - -v/ — sm 



R 



(3) 



where R is the distance to the source, is the viewing angle of the source, and AS^ i s the 



quadrupole harmonic amplitude consisting of the following three parts 11031! 15411 15811 174tl . 

^20 - ^20 (hyd) ^ ^20 (mag) 20 (grav) ' ^ ^ 

On the right hand side, the first term is related to anisotropic kinetic energies which we refer to 
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as hydrodynamic part, 
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where is effective density p, = p + (e + p + |Z7p)/c^ with p,e,p and fe^ = bjb' representing 
the baryon density, internal energy, pressure, and the magnetic energy, respectively, p = cos 6 
is the direction cosine and W - 1/ ^/T^-v'H^ is the Lore ntz boost factor with v, representing the 
spacial velocities in the spherical coordinates (see lll74ll for more details). The second term is 
related to anisotropy in gravitational potentials which we call as the gravitational part. 
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where <t) is the gravitational potential of self-gravity, and finally the third term is related to 
anisotropy in magnetic energies that we refer to as magnetic part. 
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The middle right panel of Figure [14] shows contributions to the total GW amplitudes (equa- 
tion (6)) for the strongly magnetized model (corresponding to the top left panel in Figure fT4l) . 
in which the left-hand-side panels are for the sum of the hydrodynamic and gravitational part, 
namely log (± [/fo^^^y^j + fio(sr^y)]) (1^^ top(-i-)/bottom(-)(equations and the right-hand- 

side panels are for the magnetic part, namely log (± ffo^m^g)) (right top(-i-)/bottom(-)) (e.g., equa- 
tion (fTOl)). By comparing to the middle left panels, it can be seen that the positive contribution is 
overlapped with the regions where the MHD outflows exist. The major positive contribution is 
from the kinetic term of the MHD outflows with large radial velocities (e.g., +p,W'\'r^ in equa- 
tion (|6]l). The magnetic part also contributes to the positive trend (see top right-half in the middle 
right panel (labeled by mag(-(-))). This comes from the toroidal magnetic fields (e.g., +^0^ in 
equation (fTOb). which dominantly contributes to drive MHD explosions. 

The bottom panels in Figure [14] show the GW spectra for a pair of models that does or does 
not have the increasing trend (left and right) as in the top panels. Regardless of the increasing 
trend, the peak amplitudes in the spectra are rather broad-band around 100 - 1000 Hz. On the 
other hand, the spectra for lower frequency domains (below ~ 100 Hz) are much larger for 
the model with the increasing trend (left) than without (right). This reflects a slower temporal 
variation of the secular drift inherent to the increase-type waveforms (e.g., top panels in Figure 
\T4i . Similar to the neutrino GWs mentioned in section [TT] it is true that the GWs in the low 
frequency domains are difficult to detect due to seismic noises, but a recendy proposed future 
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Figure 14: Top panels show examples of gravitational waveforms with the quasi-monoto nicall y increasing trend (left) or 
not (right panel) obtained in the 2D SRMHD simulations in Takiwaki and Kotake (201 1) |174|| . Initial rotation parameter 
is both set to be /J = 0.1 %, while the precollapse magnetic field is taken as lO'^ G (left) and lO" G (right), respectively. 
The total wave amplitudes are shown by the red line, while the contribution from the magnetic fields and from the 
sum of hydrodynamic and gravitational parts are shown by blue and green lines, respectively. The middle panels show 
various quantities for analyzing the increasing trend taking the top-left-panel model as a reference. Left panel shows the 
distributions of entropy [/.g/baryon] (left-half) and logarithm of plasma fS (right-half) at lOOms after bounce. The white 
arrows (left-half) show the velocity fields, which are normalized by the scale in the middle left edge (0.5c). The middle 
right panel shows the sum of the hydrodynamic and gravitational parts (indicated by "hyd and grav" in the left-hand side) 
and the magnetic part (indicated by "mag" in the light-hantl side), respectively. The top and bottom panels represent the 
positive and negative contribution (indicated by (+) or (-)) to A^^ , respectively (see text for more details). The side length 
of the middle panels is 4000(km)x8000(km). The bottom panels indicate the detectability of the GW spectra for a model 
with (left) or without (right) the increasing trend, respectively. The supernova is assumed to be located at the distance of 
10 kpc. 



space int erferomet ers like Fabry-Perot type DECIGO is designed to be sensitive in the frequency 
regimes (e.g., the black line in the bottom panels). These low-frequency signals, if 



observed, could be one important messenger of the increase-type waveforms that are likely to be 
associated with MHD explosions exceeding 10^' erg. 

Concerning the microphysical aspects, the MHD simulations presented in this section are 
far behind the most up-to-date simulations in which the spectral neutrino transport is solved 
(e.g., section 12] and Table 1). This means that the GW predictions in the MHD explosions are 
also still in their infancy. One may easily guess that it is numerically challenging to solve the 
neutrino transport in the highly non-spherical environments that are inherent to MHD explosions. 
The PNS is deformed to be like a dumbbell and the regions outside the PNS become highly 
anisotropic such as by the formation of the polar funnel region after the bipolar flows pass b\l~^l. 
It is an urgent task to sophisticate the GW predictions in the MHD explosions in the same level 
as those in the neutrino-driven explosion models. 

3.2.2. Non-axisymmetric instabilities 



As the textbook by Chandrasekhar (1969) 1117711 says (see also 11781 l24ll ). rapidly rotating 



compact objects can be subject to non-axisymmetric rotational instabilities when the ratio of 
rotational to gravitational potential energy (r/|W| = /I) exceeds a certain critical value. Since 
the growing instabilities carry the object's spheroidal into a triaxial configuration with a time- 
dependent quadrupole(or higher) moment, strong GW emission can be expected (see Fryer and 



New (201 1) 112311 for a detailed review of the uncovered topics in this article). 

Probably the best understood type of instability is the classical dynamical bar-mode instabil- 
ity with a threshold of /3bai ^ 0.27. Rampp et al. (1998) Ill79i1 reported the first 3D hydrodynamic 
core-collapse simulations to study the growth of the bar mode instability and their impact on the 
GW emission. The initial condition for their study was based on the configuration at several mil- 
liseconds before core bounce in the rapidly rotating 2D models of Zwerger et al. ( 1 996) il80ll . In 
addition to the configuration, they imposed low mode (m - 3) density perturbation and followed 
the growth of the instability, where m stands for the azimuthal quantum number They observed 
the three clumps merged into a bar-like structure due to the growth of the non-axisymmetric in- 
stability. In fact, their models are rapid rotators whose values of T/\W\ exceeds the critical value 
beyond ySbar, beyond which MacLaurin spheroids become dynamically unstable again triaxial 
perturbations. However, they found that the maximum GW amplitudes were only ~ 2% different 
from the 2D cases by Zwerger et al 118(1]. 3D SPH simulations of a rapidly ro tating model with 



the initial value of r/|W|init of ~ 3% by Fryer and his collaborators 181 1 reached a similar 



conclusion. Full GR 3D simulations by Shibata and his collaborators, on the other hand, pointed 
out that the maximum GW amplitudes for models with rapid rotation (1 < r/|W|init ^ 3%) can 
be en hanced by a factor of 10 than the ones predicted in the Newtonian studies mentioned above 



1 18211 . They discussed that the enhancement of the self-gravity due to the GR effects results in a 
more efficient spin-up of the core thus leading to the growth of the non-axisymmetric instability. 
The key question is whether such a precollapse rapid rotation (r/|W|init > 1%) can be realized 
or not. 

More recently, Ott et al. (2007) iflilfisl and Dimmelmeier et al. (2008) llTToll reported a 



^^In such a case, the ray-b y-ray approximation would be not good. Instead, the MGFLD scheme or a recently proposed 
Ml closure technique Il66ll would be much better. 
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systematic study in which a wide variety of the state-of-the-art stellar evolution model^ was 
employed in their 3D full GR and 2D CFC core-collapse simulations. Their results presented a 
strong evidence that the postbounce core, even in their extreme models, do not reach values of yS 
close to ySdyn during collapse, bounce and during early postbounce time. The PNS rotation could 
reach ySdyn if it would keep gravitationally contracting with conserving its angular momentum in 
their cooling phase (lasting on a timescale of ~ minutes). However in this case, the threshold of 
the secular bar-mode instability (ySsec ~ 0.14) may be first satisfied (e.g., lll84ll '). For more in- 
depth reviews about this topic with their potenti al GW emission mechanisms, see recent reviews 
of Andersson (2003) ^ and Kokkotas (2008) iflSS I. 
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Figur e 15: Time evolution of various quantities in a 3D model with r/|W|j„it = 0.2% calculated by Ott et al. (2005) 
llSdl . Time is measured from the epoch of core bounce ti,. Top panel shows that the amplitude of the m = 1 mode 
precedes that of the m = 2. Middle panel shows the time evolution of T/\W\ and the core's maximum density. It is shown 
from the panel that after the epoch of t - tj, ~ 100 ms, when the m = 2 mode begins to be ampHfied, the transfer of the 
angular momentum becomes active which results in the increase of the maximum density and the decrease of the T/l W|. 
The bottom panel shows the gravitational strain at the distance to the source r as viewed down the rotational axis (solid 
curve) and as viewed along the equatorial plane (dotted curve). One can see that the waveform traces the time evolution 
of the m = 2 mode. Note in the panel that rh =100 cm con'esponds to /i ~ 3 X 10"~' for a galactic supernova. This plot 
is by courtesy of Ott and the coauthors. 



In contrast to the high r/|W| instabilities mention ed above, recent work, so me of which has 
been carried out in idealized setups and assumptions 19111 and later also in 

more self-consistent core-collapse simulations 1 186, 1 1 IT 1 12ll . suggest that a diff'erentially rotat- 
ing PNS can become dynamically unstable at much lower T l\W\ as low as < 0.1. Despite clear 
numerical evidence for their existence, the physical origin of the low-r/|W| instability remains 
unclear. However it has been suggested [189, 192] that the instabilities are associated with the 
existence of corotation points (where the pattern speed of the unstable modes matches the local 
angular velocity) inside the star and are thus likely to be a subclass of shear instabilities 



^^also by changing a huge parameter space spanned by initial rotation rate, degree of differential rotation, and different 
EOSs. 

•^'Note that corotation resonance has been long known to the key ingredients in the accretion disk system, such as the 



Papaloizou-Pringle instability [19: 
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Figure [TS] shows various quantities for a non-magnetized 20 Mq model of Ott et al. (2005) 
that experiences the low-r/|W| instability. They investigated the gro wth o f the non-axisymmetric 
structure until the rather later phases (> 100 msec) after bounce |l86ll. They foun d that the 
growth of the m - 1 mode, the so-called one- armed instability il94ll87l [l88 l [l92ll. precedes 
the growth of the bar-mode (m = 2) instability (see the top panel in Figure [TSl). where m denotes 
the azimuthal quantum number. They pointed out that the initial rotation rate can be as small 
as r/|W|init = 0.2% for the sufficient gravitational radiation enough to be detected by the future 
detectors. This value is indeed much smaller than the one previously assumed for igniting the 
growth of the classical bar-mode instability. 



More recently, Scheidegger et al. (2008,2010) 111 1 U ll 1211 investigated this issue by more 



realistic 3D MHD models that followed consistently the dynamics from core-collapse, through 
bounce, up to ~ 200 ms postbounce in which a treatment for neutrino transport based on a partial 
implementation of the IDS A scheme 1 9^ 3 was incorporated. Top panels in Figure [16] show the 
vorticity distribution in the equatorial plane selected from their 3D MHD models, in which the 
major difference originates from the inclusion of the neutrino transport effect (right) or not (left). 
Regardless of the model difference, a two-armed (m - 2) pattern is clearly seen (at ~ 10 ms after 
bounce for the plot) and then the spiral flows develop for more than several hundreds millisec- 
onds later on. The middle panels display the quadrupole amplitude of the A+ and Ax along the 
pole for their 3D MHD model with (right) or without the neutrino cooling (left). The most strik- 
ing feature is that the model with neutrino cooling (right) shows 5 ~ lOx bigger maximum GW 
amplitudes due to the nonaxisymmetric dynamics compared to their counterparts that neglect 
neutrino cooling (left). They pointed out that the neutrino cooling in the postbounce phase leads 
to a more condensed PNS with a shorter dynamical timescale and also with much more matter 
enclosed in the unstable region (see the compactness of the spiral flows in the top right com- 
pared to the top left panel in FigurefTSI), leading to much more powerful GW. Another important 
message from the middle panels is that the gravitational waveforms from the nonaxisymmetric 
dynamics generally show narrow-band and highly quasi-periodic signals which persist until the 
end of simulations. For the periodic signals, the effectively measured GW amplitude scales with 
the number of GW cycles as /leff <x h ^/N. This is the reason why the peak amplitude in the 
GW spectrum (the bottom panel in Figure [T6]l exceed 10"^° for the galactic source, making the 
chance of detection quite higher. In fact, the wave amplitude from the short-duration bursts near 
bounce is generally smaller than 10"^° for the galactic source (e.g., the amplitude near 1 kHz in 
the bottom panels in Figure 17). 

For detecting (potentially) the most powerful GW signals, it is therefore of crucial importance 
to understand the properties of the non-axisymmetic instabilities. However the numerical diffi- 
culty to follow a long-term postbounce evolution in 3D is the main hindrance at present. While 
the angular momentum is continuously brought in the iron core from stellar envelope that rotates 
with higher angular momentum, the shear energy in the vicinity of the corotation points will 
continue to be redistributed or even dissipated by some mechanisms including the MRI. In the 
mentioned 3D simulations, neutrino heating is not taken into account. Neutrino-driven convec- 
tion would affect the growth of the non-axisymmetric instabilities as well as through the spiral 
SASI modes. The general relativity that has been mostly treated in a post-Newtonian manner 
should affect the growth rate. At first, we could start to address this question by performing para- 
metric explosion models in 3D MHD (e.g., sections |3. 1 . Il andl 3. 1 .2l i. But for a more quantitative 



'see section|2]for more details of the scheme. 
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Figure 16: Top panels depict the vorticity distributi on in the equatorial plane of the 3D MHD models of R4STAC (left) 
and R4E1FCl (right) in Scheidegger et al. (2010) fTTll (at ~ 10 ms after bounce). A two-armed m = 2 pattern behind 
the stalled shock is clearly seen (the side length of the plot is (300km)-) either for their 3D model without (top left) or 
with the deleptonization effect (top right). The middle panels display the quadrupole amplitude of the A+ and Ax along 
the pole for a given 3D MHD model without (left) or with the deleptonization effect (right). The bottom panel shows 
spectral energy distributions of the GW signal emitted along the polar axis for their representative 3D model experiencing 
the low-T/l W| instability (model R4STCA) at a distance of 10 kpc. These plots are by courtesy of Scheidegger and the 
coauthors. 35 



discussion, the full 3D GRMHD simulation with an appropriate neutrino transport is needed also 
in this case, towards which the supernova modelers will keep their efforts (as they have done 
since time immemorial) and more recently people in the numerical relativity are also joining in 
the efforts. 



4. Summary and Concluding remarks 
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Table 3: Illustrative summarizing the relation of the potential explosion mechanisms (horizontal column) and their emis- 
sion processes (columns colored by yellow) and GW signatures (columns colored by red) obtained so far in "2D" and 
"3D" simulations (vertical direction). In each column, the waveform (left) and the GW spectrum (right) are shown for 
some representative models, in which top left (2D, neutrino mechanism), top right (2D, MHD mechanism), bottom left 
(3D, neutrino mechanism), and bottom right (3D, MHD mechanism) are again taken from Figures I4l ll4| [8l and ll6l (see 
each figure and their ca ption for more details). Figures in th e top left and bottom right columns are taken from the results 
by Miiller et al. (2004) fl^ and Scheidegger et al. (2010) fTTal by their courtesy. 



The aim of writing this article was to provide an overview of what we currently know about 
the possible GW signatures of core-collapse SNe predicted by a number of extensive numerical 
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simulations especially in the last ~15 years. In addition to the long-studied bounce GW signals, 
multiple new physical inputs have been proposed mostly in the last decade. Among them, we 
focused on the GW signatures produced by prompt convection, non-radial flows inside the PNS 
and in the neutrino-driven bubble, the activity of the S ASI, and by asymmetric neutrino emission 
in the context of the neutrino-driven explosion models in section IXTl and by asymmetries asso- 
ciated with the effects of magnetic fields and non-axisymmetric rotational instabilities in section 
13.21 As mentioned in section |2l our understanding of the physics to unravel these waveform 
features has been progressing in accordance with the sophistication of the numerical simulations 
on which SN modelers have been putting a huge eflrort for long. 

We wish that from Table 3, one might see a very rough illustrative summary of this article that 
describes a possible association between the potential explosion mechanisms (horizontal) and the 
emission processes with their GW signatures (columns colored by yellow and red) obtained in 
the recent multi-D simulations (vertical). At one sight it is a very good news for us that there 
is a clear correspondence between the explosion mechanisms and their GW signatures. This 
means that the GW observation could provide an important probe into the potential explosion 
mechanisms (see also Ott (2009) lH^lO Although the criterion that bifurcates between the 
neutrino-driven mechanism and the MHD mechanism is still gray, rotation and magnetic fields 
in the precollapse cores should hold the key importance. 

In the neutrino heating mechanism, the whole story may be summarized as follows. If a 
precollapse iron core has a "canonical" rotation rate as predicted by stellar evolution calculations 
(f2o 0.1 rad/s), the collapse-dynamics before bounce proceeds spherically and the structures 
interior to the PNS are essentially spherical at this epoch. The stalling bounce shock gives rise 
to the GW signals of the prompt convection (see the inset of the top left panel in Table 3, and 
also Figure 1 in section [TTI for more details). Typically later than ~ 100 ms after bounce, con- 
vective overturns as well as the SASI in the postshock region become much more vigorous with 
time. Since anisotropics of the neutrino flux as well as matter motions in this epoch are gov- 
erned by the non-linear hydrodynamics, the GWs signatures change also stochastically with time 
(sections 13.1.11 and 13 . 1 .2l i. This feature was already known in the first-principle 2D simulations 
(indicated by "stochastic" in Table 3, e.g., section 13.1.3 1 for more details). Furthermore recent 
exploratory 3D simulations revealed that the stochastic nature becomes much more strong be- 
cause the explosion geometry changes in all directions (section [3.1.4l l. Note that the stochastic 
nature could become weaker in the presence of rotation as a result of the spiral SASI (section 
I3.1.4I I. however the first-principle 3D simulations are apparently needed to draw a robust conclu- 
sion. Contributed by neutrino GWs dominant in the lower- frequency bands (< 100 Hz), the total 
GW spectrum tends not to be a single peak but rather flat over a broad frequency up to ~ IkHz. 
For detecting these GW signals for a Galactic supernova source with a good signal-to-noise ra- 
tio, we need next-generation detectors such as the advanced LIGO, LCGT, and the Fabry-Perot 
type DECIGO (e.g., the left panels in Table 3 (right-hand)). As a side remark, neutrino signals 
at this non-linear phase change also stochastically with time, which could be detectab le for a 
Galactic supernova in a currently running detector such as by IceCube l87l Il28l 119511 . When 
material behind the stalled shock successfully absorbs enough neutrino energy to be gravitation- 
ally unbound from the iron core, a powerful neutrino-driven explosion with its explosion energy 
exceeding 10^' erg will be obtained (section|2]l. In this case, the GW waveforms would imprint 



^'it should be also remembered that galactic core-collapse SN event is a quite rare event happening at a rate of one in 
a few decades. As theoreticians, we sincerely hope to be able to make clear the supernova mechanism and the related 
issues before observations tell us ! 
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the information of the blast morphologies (e.g., section U. l.ll i. 

If a precollapse core rotates enough rapidly (typically initial rotation period less than ~ 4 s) 
with strong magnetic fields (higher than ~ 10" G) imposed initially, the MHD mechanism can 
produce bipolar explosions along the rotational axis predominantly by the field wrapping mech- 
anism (section l3T2ll^^[ The GW signals in energetic MHD explosions exceeding 10^' erg, are 
characterized by a burst-Uke bounce signal plus a secularly growing tail (indicated by "Bounce" 
with "tail" signals in the top right panel in Table 3, also section I3.2l i. Likewise the neutrino 
GWs in the case of neutrino-driven explosions, the tail component makes the total GW spectrum 
broad-band. Going beyond 2D, recent 3D MHD simulations show a strong evidence that nonax- 
isymmetric rotational instabilities can be a source of strong GW emission due to its periodicity 
and bar-mode character (see a sharp peak near ~ IkHz in the bottom right panel in Table 3). 
Note that the pole to equator anisotropy of the shock propagation in the MHD explosions could 
lead to a sudden decrease in the SuperKamiokande events through the MSW flavor conversion 
(~ 2500 events for a Galactic source). A planned joint analysis of neutrino and GW data il96ll 
would provide a powerful probe especially into the MHD-driven explosions that imprint these 
peculiar signatures distinct from other candidate explosion mechanisms. Note also that MHD 
explosions could be a possible r-process cite because the mass ejection from the iron core takes 
place in much shorter timescales than in the delayed neutrino-heati ng rnecha nism, which pre- 
vents the ejecta from being proton-rich by neutrino capture reactions 0197 , [l98lF[ As mentioned 



in the final part of section 13.21 these interesting possibilities have been proposed so far by the 
numerical simulations with a crude treatment of neutrino transport, so that they are now awaited 
to be re-examined. 

One important notice here is that the explosion energies obtained in some of the first-principle 
2D simulations (Table 1 in section |2]i are underpowered by one or two orders of magnitudes to 
explain the canonical supernova kin etic e nergy (~ 10^' erg). Moreover, the softer nuclear EOS, 
such as of the Lattimer-Swesty (LS) 120111 EOS with an incompressibility at nuclear densities, K, 
of 180 MeV, has been employed in all the 2D simulations that succeeded in producing neutrino- 
driven explosion^ On top of a striking evidence that favors a stiffer EOS based on the nuclear 
experimental data {K - 240 + 20 MeV, 120211') . the soft EOS may not account for the recently 
observed massive neutron star of ~ 2Mq 1120311 [see the maximum mass for the LSI 80 EOS in 



204l 120511 . What is then missing furthermore? To seek the answer is equal to making the GW 



prediction more accurate. We may get the answer by going to the first-principle 3D simulations 
or by ta king i nto account new ingredients, such as exotic physics in t he core of the protoneutron 
star lf76[ l206ll . viscous heating by the magnetorotational instability 1 207 , 208 1. or energy dissi- 
pation via Alfven waves ll209ll . General relativity (GR) so far treated in a very approximate way 
in multi-D models might help explosions if the advantage of the GR (i.e., higher neutrino ener- 
gies) dominates over the disadvantage (i.e., shorter advection timescale). 3D effects might help 
neutrino-driven explosions as have been recently pointed out by 3D simulations with a spectral 
neutrino transport lf79ll . Assisted by a growing computer power, all these issues will be tested 
soon by the forth-coming first-principle simulations. 



If the MRI can be sufficiently resolved in global simulations, MHD outflows might be produced even for more 
weakly magnetized cores because the MRI can exponentially amplify the initial magnetic fields to a dynamically impor- 
tant strength within several rotational perio ds (e.g., collective references in Obergaulinger et al. (2009) llVlll ). 

See collective references in IT99II2OOI1 for other plausible r-process cites. 
^''2D explosion was obtained in the MPA simulations for a non-rotating 1 1.2 M© for the Shen EOS (K = 281MeV) 
(H-T. Janka, private communication). 
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In addition to these advances in numerical simulations, the physical understanding of funda- 
mental problems in the supernova theory is progressing at the same pace. What determines the 
saturation levels of S ASI? A careful analysis on the parasitic instabilities has been reported to 
answer this question 1210( 1. What determines mode couplings between the sm all-s cale convec- 
tion eddies and the large scale SASI modes? To apply the theory of turbulence 121111 might put a 
milestone to address this question. The physical understanding of these phenomena should also 
give us a deep understanding of the GW features from the SASI and convection. Although it is 
also numerically challenging to resolve the small scale MRI modes in global simulations, it is 
very interesting to see their outcome in the neutrino-driven mechanism with rapid rotation. Be- 
sides the two representative EOSs of LS and Shen, new sets of EOSs have been recently reported 
ll95il96[ 971. Can we extract the information of the different EOSs by the GW observation ? All 



of these questions may not be easy to answer immediately, but they should indeed provide us a 
precious chance to solve them if not for the GW astronomy. 

If the neutrino-heating mechanism (or some other mechanisms) fails to explode massive stars, 
central PNSs collapse to BHs. Recent GR simulations ll212ll indicate that the significant GW 
emission is associated at the moment of the BH formation, which can be a promising target of 
the advanced LIGO for a Galactic source. As pointed out in 112130 . disappearing neutrino signals 
could also tell us the epoch of BH formations for a Galactic source. When quarks and pions could 
appear in a collapsing core, one may observe a neutrino burst induced by the recollapse of the 
PNS, which is triggered by the sudden softening of the EOS i76il . The intervals betw een bounc e 
and the BH formation depend on the details of exotic physics in the super-dense core I2l4l215ll . 
All of these observational signatures co uld b e an important probe into the so-called dense QCD 
region in the QCD phase diagram (e.g., 1121611 for recent review) to which lattice calculations are 
hardly accessible at present. 

Finally, what is about the story if the MHD mechanism fails ? The central cores with signifi- 
cant angular momentum would c ollapse int o a BH. Neutrinos emitted from the accretion disk heat 
matter in the polar funnel region 12171121811 or strong magnetic fields in the cores of order of 10'^ 
G play also an active role both for drivi ng the rnagnet o-driven jets and for extracting a significant 



amount of energy from the BH (e.g.. 11331 
often referred as collapsar fl22U l222l |223| 



and see references therein). This picture. 



has been the working hypo thesis as a central engine 

of long-duration gamma-ray bursts for these 20 years (see references in ll224il219[|220l.l22ll226ll 
for other candidate mechanism including magnetar models). It is indeed a big issue whether the 
formation of outflows predominantl y proceeds via the neutrino-heating mechanism or MHD pro- 
cesses (see collective references in 12271 12281 19211 for example). Comparing to various findings 
illustrated in Table 1, much little things are known about the BH-forming supernovae and about 
the collapsar This may be partly because the requirement for the numerical modeling to this 
end is highly computationally expensive, which necessitates not only GRMHD simulations for 
handling the BH formation, but also multi-angle neutrino transfer for treating highly anisotropic 
neutrino radiation from the accretion disks. 

To get a unified picture of massive stellar death in view of GW astronomy, we need to be 
able to draw a schematic picture (like Table 3) also in the case of the BH forming supernovae. A 
documentary film recording our endeavours to make our dream of "GW astronomy of massive 
stellar core-collapse" come true, seems not to come to an end immediately and is becoming 
even longer by taking into account the BH forming supernovae. In addition, a new horizon is 
now being opening up towards a multi-messenger astronomy by which we could have much 
more information to decipher the central engines of massive stars from the combined analysis 
between GWs, neutrinos, and electromagnetic messengers. For example, optical observations 
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using a large-aperture telescope such as the Subaru telescope with a detailed spectropolarimetric 
technique (e.g., lf28ll ) is expected to give us an indication of the explosion geometry, which 
could be also helpful to constrain the asphericity of the central engines. Needless to say, it is 
of crucial importance to ac curately d etermine nucleosynthesis in the SN ejecta for modeling 



supernova light curves (see, J229ll23(ill for recent reviews). Note that these studies are now able 



to be updated by using recent results obtained in the first-principle simulations as mentioned in 
section|2] Here we would like to lay down our pen hoping that in the near future a great progress 
will be made in our understanding of a number of exciting and unsettled issues raised in this 
article. 
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